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Abstract. Let V = (x^^pij{x), i — 1, . . . ,n, j = 1, . . . ,Ni) be a space of quasi- 
polynomials in x of dimension N =^ Ni + ■ ■ ■ + Nn ■ The regularized fundamental 
differential operator of V is the polynomial differential operator "^^^^^ Aj\i^i{x){xdxy 
annihilating V and such that its leading coefficient Aq is a monic polynomial of the 
minimal possible degree. Let U — (z1^ <lab{u), a — 1, ... ,m, b — 1, . . . , Ma) be a space 
of quasi-exponentials in u of dimension M = Mi + - ■ ■+Mm. The regularized fundamen- 
tal difference operator of U is the polynomial difference operator X]f=o ^A^-j(^) ('''")' 
annihilating U and such that its leading coefficient Bq is a monic polynomial of the 
minimal possible degree. Here (ru/)(u) = f{u + 1). 

Having a space V of quasi-polynomials with the regularized fundamental differential 
operator D, we construct a space of quasi-exponentials U = ( z" qab{u) ) whose regular- 
ized fundamental difference operator is the difference operator M*Ajv-i(T„). The 
space U is constructed from by a suitable integral transform. Similarly, having U 
we can recover by a suitable integral transform. Our integral transforms are analogs 
of the bispectral involution on the space of rational solutions to the KP hierarchy |^ . 

As a corollary of the properties of the integral transforms we obtain a correspon- 
dence between solutions to the Bethe ansatz equations of two (flfA?, sCm) dual quantum 
integrable models: one is the special trigonometric Gaudin model and the other is the 
special XXX model. 

1. Introduction 

Let V = { x^^pij{x), i = 1, . . . ,n, j = 1, . . . , Ni) he a. space of quasi-polynomials in x 
of dimension = A^^i -|- ■ ■ ■ -|- A^„. The regularized fundamental differential operator of 
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V is the polynomial differential operator Yli=o ^N-i{x){xdxy annihilating V and such 
that its leading coefficient Aq is a monic polynomial of the minimal possible degree. 

Let U = {z"^ (lab{u), a = 1, . . . ,m, 6 = 1,..., Ma) be a space of quasi-exponentials in 
u of dimension M = Mi + ■ ■ ■ + Mm- The regularized fundamental difference operator of 
U is the polynomial difference operator YlfLo BM-i{u){TuY annihilating U and such that 
its leading coefficient Bq is a monic polynomial of the minimal possible degree. Here 
{rj){u) = f{u + l). 

Having a space V of quasi-polynomials with the regularized fundamental differential 
operator AN-i{x) {xdx)\ we construct a space of quasi-exponentials U = {z^ qab{u) ) 
whose regularized fundamental difference operator is the difference operator 
MMAr_i(r„). The space U is constructed from by a suitable integral transform. 

Having a space U = {z^ qahiu) ) of quasi-exponentials with the regularized fundamen- 
tal difference operator -^A^-i('w)(T„)*, we construct a space of quasi-polynomials 

V = {x^^Pij{x)) whose regularized fundamental differential operator is the differential 
operator YldLo^^^M-ii^dx)- The space V is constructed from f/ by a suitable integral 
transform. 

Our integral transforms are analogs of the bispectral involution on the space of rational 
solutions to the KP hierarchy jWj. 

As a corollary of the properties of our integral transforms we obtain a correspondence 
between solutions to the Bethe ansatz equations of two {qIn^Q^m) dual quantum in- 
tegrable models: one is the special trigonometric Gaudin model and the other is the 
special XXX model. 



Example. Let n = (ni,n2) and m = (mi, 7712) be two vectors of nonnegative integers 
such that ni + n2 = mi + m2. Let d be the number of integers i such that max (0, n2 — 
mi) < i < min (m2, n2). Let z = {zi, Z2), A = (Ai, A2) G be two points with distinct 
coordinates. 

Consider two systems of algebraic equations. The first system is the system of equa- 
tions 

, . Ai — A2 — 1 X — ^ \ — ^ 2 

(LI) -^—^ — + Er^ = 0' ^ = i^---^^2, 

a=l j=l, jj^i J 

with respect to the unknown numbers ti, . . . , t^j. The system is symmetric with respect 
to the group of permutations of ti, . . . , t„2. One can show that for generic z, A the 
number of S^j-orbits of solutions of system (jl.ll) is equal to d. This system is called the 
system of the Gaudin Bethe ansatz equations, see Section 110.51 
The second systems of equations is the system 

2 . r?i2 

(1-2) II 7 II TT = ~ ' a = l,...,m2, 

1=1 b, b^a 
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with respect to the unknown numbers Si, . . . , Sm^. The system is symmetric with respect 
to the group of permutations of Si, . . . , Sm2- One can show that for generic z, A the 
number of -orbits of solutions of system ()1.2|) is equal to d. This system is called 
the system of the XXX Bethe ansatz equations, see Section 110.51 

Problem. Establish a correspondence between orbits of solutions to systems fll.ip and 

(El. 

We give two solutions to this problem. 

The first solution. To system we assign the vector space (L^j L^a) ['^i, '"'2] and 

four commuting linear operators acting on this space. Here (L^.^ ^ Lm2)[ni,n2] denotes 
the weight subspace of weight [ni,n2] of the tensor product of g[2"^odules with highest 
weights (mi,0) and (m2,0), respectively. The space (i^mi ® 

) [ni, n2] is of dimension 

d. The linear operators are denoted by 

Hf{Xi + ni, A2 + n2, z) , Hf{Xi + rii, A2 + na, z) , 

Gf (Ai + A2 + n2, z) , Gf{\i + rii, A2 + n2, z) 

and called the Gaudin KZ and dynamical Hamiltonians. To each orbit of solutions 
of system p.ip . the Bethe ansatz method assigns an eigenvector of the commuting 
Hamiltonians. The constructed Bethe vectors form a basis of this vector space. 

To system ()1.2|) . we assign the vector space (L„^ (8)L„2) [?tii, 777.2] of the same dimension 
d and four commuting linear operators acting on this space. The linear operators are 
denoted by 

Hf{z, Ai + ni, A2 + 772) , H^{z, Ai + rzi, A2 + 772) , 

Gf(z, Ai + 77i, A2 + 7^2) , (2, Ai + rii, A2 + 772) 

and called the XXX KZ and dynamical Hamiltonians. To each orbit of solutions of 
system fjl.2|) . the Bethe ansatz method assigns an eigenvector of the commuting Hamil- 
tonians. The constructed Bethe vectors form a basis of this vector space. 

There is a natural isomorphism of the vector spaces {Lm^ ® ) [^1 , ''^2] and ® 
L^j) [77ii, 77^2], which identifies the Hamiltonians: 

iff (Ai + 771, A2 + 772, Z) = (Z, Ai + 771, A2 + ^2) , 

Gf (Ai + 77i, A2 + 772, Z) = Hf{z, Ai + 77i, A2 + 772) , 

for i,a = 1,2. This isomorphism is called the (0I2, 0[2)-duality, see |TV5j . Under the 
duality isomorphism the Bethe vectors are identified and this identification gives a cor- 
respondence between the orbits of solutions of systems ()1.1|) and p.2|) . 

The second solution. To the orbit of a solution {ti, . . . ,tn2) of system ()1.1|) . we assign 
the polynomial p2{x) = 11^=1 ~ ^i) differential operator 

D = x\x - z,){x - Z2){d. - ln\- ^^"=^^ °^ )) (9. - ln'(x^^p2)) , 

P2 
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where dx = d/dx and for any function /, In' / denotes the logarithmic derivative /'//. 
Clearly the differential equation Df{x) = has a solution x'^'^ p2{x). 
We show that D can be written in the form 

D = Ao{x){xdx)^ + Ai{x)xdx + A2{x) 

where Aq, Ai, A2 are polynomials in x of degree not greater than two. Then we consider 
the second order difference equation 

(mMoK) +uAi{t^) +A2{r^))g{u) = 

for the unknown function g{u). It turns out that this difference equation has a solution 
of the form Z2q2{u), where q2{u) = Y[T=i ('^ ~ ^a), and the roots Si, . . . , 5^,2 form a 
solution of system ()1.2|) . This construction gives a second correspondence between the 
orbits of solutions of systems and ()1.2j) . 

We show that the two described constructions (the first, based on the (gig, 3l2)"duality, 
and the second, which uses the differential and difference operators) give the same cor- 
respondence between solutions of systems and ()1.2|1 . 



This paper is a development of results of the paper |MTVlj , in which we presented an 
integral transform giving an involution on the space of quasi-exponentials, the involution 
which corresponds to the bispectral involution of G.Wilson in |Wj . 



The paper has the following structure. In Section El we discuss spaces of quasi- 
polynomials and their fundamental differential operators. In Sectional we discuss spaces 
of quasi-exponentials and their fundamental difference operators. In Section |3J we de- 
fine integral transforms establishing a duality between spaces of quasi-polynomials and 
quasi-exponentials. In Sectional we introduce special spaces of quasi-polynomials and 
quasi-exponentials. We introduce special integral transforms relating the special spaces. 
In Section we discuss relations between special spaces and solutions of the Bethe 
ansatz equations. The special spaces of quasi-polynomials correspond to solutions of 
the Bethe ansatz equations in (trigonometric) Gaudin models. The special spaces of 
quasi-exponentials correspond to solutions of the Bethe ansatz equations in XXX mod- 
els. The number of solutions of the Bethe ansatz equations is discussed in Section |71 
We describe the Gaudin and XXX models in Section |H1 The {QlN^dhi) duality between 
the (trigonometric) Gaudin and XXX models is defined in Section 121 In Section El we 
formulate a conjecture about the bispectral correspondence of Bethe vectors under the 
i&^N^dhi) duality. In Section (TUl we prove the conjecture for = M = 2. In Section 
im we discuss the relation of our integral transforms to the bispectral correspondence 
of suitable Baker- Akhieser functions for differential and difference operators. 

Authors thank Y. Berest and O. Chalykh for discussions. 
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2. Spaces of quasi-polynomials 

2.1. Definition. Let p G C[x] be a polynomial, A a complex number. The function 
x^p is called a quasi-polynomial in x. The quasi-polynomial is a multi-valued function. 
Different local uni-valued branches of the function differ by a nonzero constant factor. 

Let A^i, . . . , Nn be natural numbers. Set N = Ni + ■ ■ ■ + iY„. For i = 1, . . . , n, let 
< nji < ■ ■ ■ < TijjVj be a sequence of positive integers. For i = 1, . . . ,n, j = 1, . . . ,Ni, 
let Pij e C[x] be a polynomial of degree Uij. 

Let Ai, . . . , A„ e C be distinct numbers such that Aj — Aj ^ Z for i ^ j. 

Denote by V the complex vector space spanned by functions x^'-pij{x), i — 1, . . . ,n, 
j = 1. . . . , Ni. The dimension of V is N. 

The space V is called the space of quasi-polynomials. 

We say that the space V is non-degenerate if for any i = 1, . . . ,n and any j = 1, . . . , Ni, 

• there exists a linear combination of polynomials Pii,Pi2, ■ ■ ■ ,PiNi which has a root 
at X = of multiplicity j — I, 

• the space V does not contain the function x^^'^'^'^ . 

2.2. Exponents. Let V be the space of quasi-polynomials. 

For z E C*, define the sequence of exponents of V at z as the unique sequence of 
integers, e = {ei < ■ ■ ■ < ejy}, with the property: for i = 1, . . . , N, there exists f & V 
such that / has a root a.t x = z of multiplicity e,. 

We say that 2; e C* is a singular point of V if the set of exponents of F at 2; differs 
from the set {0, . . . , — 1}. The space V has at most finitely many singular points. 

Let {zi, . . . , Zm) be the subset of C* of all singular points of V . For a = 1, . . . , m, let 

{Q<---<N -Ma-KN -Ma + mal<---<N -Ma + TUaMa } 

be the exponents of V aX Za. Here < mai < • • • < rUaMa ^iid is an integer such 

that 1 < Ma < A^. Set Af = Ml H h M„,. 

We say that ^ is a space of quasi-polynomials with data 

Dy = {n,Ni,nij,\i, m,Ma,mab,Za} 
where i — 1, . . . ,n, j — 1, . . . , Ni, a — 1, ... ,111, b — 1, ... , Ma. 

2.3. Fundamental differential operator. For functions /i, . . . , /j of one variable, de- 
note by Wr (/i, . . . , /j) their Wronskian, that is, the determinant of the i x i-matrix 
whose j-th row is fj, fj^\ . . . , /j* ^\ 

Define the Wronskian of V, denoted by Wry, as the Wronskian of a basis of V. The 
Wronskian of V is determined up to multiplication by a nonzero number. 

Lemma 2.1. Let V be a non- degenerate space of quasi-polynomials, then 

m Ma n Ni 

a=l 6=1 i=l j=l 
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The lemma is proved by analyzing the order of zeroes of the Wronskian of V and its 
asymptotics at infinity. 

The monic fundamental differential operator of V is the unique monic linear differen- 
tial operator of order N annihilating V . It is denoted by Dy- We have 

Dy = + + + A = (-ir^> 

where = d/dx, Wry is the Wronskian of a basis {/i, . . . , /jv} of V, Wvy^i is the deter- 
minant of the N X iV-matrix whose j-th row is fj, fj^\ . . . , fj^ * ^\ fj^ ^'^^\ . . . , fj^^- 
For any j — 1, . . . , N, the order of the pole of Aj at x — Za, a — 1, . . . , m, does not 
exceed j. 

Lemma 2.2. Let V be a non- degenerate space of quasi-polynomials, then for a — 
1, . . . ,m, the order of the pole of Am^ o-t Ma and the order of the pole of 

Ai is not greater than Ma for i > Ma- 

The proof follows from counting orders of determinants Wy^i. 
The polynomial differential operator 



Dy = x'^Hix-Zar^ D, 



'V 

a=l 

is called the regularized fundamental differential operator of V . 

It is easy to see that the regularized fundamental differential operator Dy of V has the 
form Aq{x)x^ + Ai[x)x'^ ■ ■ ■ + Aj^{x) where Ai{x) are polynomials. Using the 
formula a;*9* = xdx{xdx — 1) ■ ■ ■ {xd^ — i -|- 1) we can present the regularized fundamental 
differential operator in the form 

Dy = Ao{x){xd,f+ A^{x){xd,f-^+ A2{x){xd,f-^+ ■■■ + Am{x) 

with polynomial coefficients i — 0, . . . , N. 

Lemma 2.3. Let V be a non- degenerate space of quasi-polynomials, then 

• WehaveAoix) = 11^=1 - -^a)"^"- 

• All coefficients Ai{x) are polynomials in x of degree not greater than M . 

• Write Dy = x^Bo{xd^) + x^-^Bi{xd^) + ■■■ + xBM-i{xd^) + BM^xd^) where 
Bo{xdx), . . . , BM{xdx) are polynomials in xd^ with constant coefficients. Then 

n Ni m n Ni 

Bo{xd,) = H Hixd, - A, - n,,) , BM{xd,) = 1[{-Za)^'' J] H^^-^- - A, - j + 1) . 

i=l j=l a=l 1=1 j=l 

• The polynomials Bq, . . . , Bm have no common factor of positive degree. 

□ 
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2.4. Conjugate Space. Let V he a space of quasi-polynomials as in Section ITTl 

The complex vector space spanned by all functions of the form Wr(/i, . . . , Jn-i) /Wry 
with fi & V has dimension A^. It is denoted by V* and called conjugate to V. 

The complex vector space spanned by all functions of the form 
f{x) x''^ YY^^^{x — Za)~^'' with / G is denoted by V'^ and called regularized con- 
jugate to V . 

Lemma 2.4. Let V he a non- degenerate space of quasi-polynomials, then 

(i) For a = 1, . . . ,m, if e are exponents of V at Za, then 

e* = {-cn -l + N < -CN-i -l + N <■■■ <-ei-l + N} 
are exponents of V* at Za and 
= {-e^v -1 + N - Ma< -e^_i - l + iV-M„<---<-ei-l + iV - M„} 

are exponents ofV'^ at Za- 

(ii) For any i = 1, . . . ,n, j = 1, . . . ,Ni, there exists f E V'^ such that the function 
x^i+j f has a nonzero limit as x ^ 0. 

(iii) For any i = 1, . . . ,n, j = 1, . . . ,Ni, there exists / G such that the function 
^Xi+riij+i-Ni j- f^^g ^ nonzero limit as x oo. 

□ 

Let D = t>e a differential operator with meromorphic coefficients. The 

operator D* = ^^ii—dxYAilx) is called formal conjugate to D. 

Lemma 2.5. Let V be a non- degenerate space of quasi-polynomials. Let Dy and Dy 
he the monic and regularized fundamental differential operators ofV, respectively. Then 
{Dv)* annihilates V* and {Dv)* annihilates V'^ . □ 

3. Spaces of quasi-exponentials 

3.1. Definition. Define the operator r„ acting on functions of u as {Tuf){u) = + 1). 

Let z be a nonzero complex number with fixed argument. Set 2;" = e"'°^. We have 
r„z" = z'^z. 

Let q G C[u] be a polynomial. The function z^q is called a quasi- exponential in u. 

Let Ml, . . . , Mm be natural numbers. Set M = Mi + ■ ■ ■ + M^. For a = 1, . . . , m, 
let < niai < ■ ■ ■ < rUaMa be a sequence of positive integers. For a = 1, . . . ,m, b = 
1, . . . , Ma, let qab G C[n] be a polynomial of degree ruab- 

Let Zi, . . . , Zm he distinct nonzero complex numbers with fixed arguments. 

Denote by U the complex vector space spanned by functions z'^qabiu), a = 1, . . . , m, 
6=1,..., Ma- The dimension of U is M. 

The space U is called the space of quasi-exponentials. 

For functions /i, . . . , /j of m, denote by Wr^*^-* (/i, . . . , /j) their discrete Wronskian which 
is the determinant of the i x -i-matrix whose j-th row is fj (u) , fj{u -\- 1) , . . . , fj{u -\- i — 1) . 
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Define the discrete Wronskian ofU, denoted by Wr^ , as the discrete Wronskian of a 
basis of U. Tlie discrete Wronskian of U is determined up to multiplication by a nonzero 
constant. 

Lemma 3.1. We have 

m 

Wr[?)(«) = S{u) n ^f'"" , 

a=l 

where S{u) is a polynomial of degree Yll=i 12j^=ii^ij + 1 ~ j)- ^ 

3.2. The frame of a space of quasi-exponentials. Let f/ be a space of quasi-expo- 
nentials like in Section ITTl Let vi, . . . ,vm be the quasi-exponentials generating U. 

For i = 1, . . . , M, let Si G <C[u] be the monic polynomial of the greatest degree such 
that the function WT^^\vjj^,Vj2, . . . ,Vj^)/Si is regular for any ji, j2, • • • , ji £ {1, • • • , M}. 

In particular, for the discrete Wronskian of U we have 

m 

Wr[?)(M) = const JJ zf"" 

a=l 

with a nonzero constant. 

Lemma 3.2. There exists a unique sequence of monic polynomials Pi{u), . . . , Pm{u) 
such that 

i i—k+1 

s^{u) = n n Pk{u+3~i) 

k=l 3=1 

fori = l,...,M. 

This lemma is an analog of Lemma 4.12 in |MV2j . 

The monic polynomials Pi{u), . . . , Pm{u) are called the frame of U . 

Proof. We construct Pi by induction on i. For i = 1, we set Pi = Si. Suppose the 
lemma is proved for alH = 1, . . . , io ~ 1- Then we set 

R{u) = n n p^M = Si,{u)/R{u). 

i=i j=i 

We have to show that Pig is a polynomial. In other words, we have to show that the 
discrete Wronskian of any io-dimensional subspace (fj^, . . . ,Vj.^) is divisible by R{u). 

Consider the Grassmannian Gr{iQ — 2, U) of {iq — 2)-dimensional spaces in U. For any 
2; G C the set of points in Gr{io — 2,U) such that the corresponding discrete Wronskian 
divided by SiQ-2 does not vanish at z, is an open set. Therefore, we have an open set of 
points in Gr{io — 2, U) such that the corresponding discrete Wronskian divided by Sig-2 
does not vanish at roots of R{u — 1). We call such subspaces acceptable. 

Therefore, we have an open set of points in Gr{io, U) such that the corresponding io- 
dimensional space contains an acceptable io — 2 dimensional subspace. Let Wi, . . . , Wi^ G 
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U be such that wi, . . . ,WiQ-2 span an acceptable space. It is enough to show that 
Wy^'^\wi, . . . , Wjg) is divisible by R{u). 

Using discrete Wronskian identities of |M V2j . we have for suitable holomorphic func- 
tions flj2,9- 



Wr('^)( 



wi, . . . ,m 



10 J 



Wr(<^)K,...,W7,„_2)(M + l) 



Sio-2{u + i)giu + 1) Si,,-2iu + i) giu + i) giu + i) 

Since the space spanned by Wi, . . . , Wj(,_2 E U is acceptable, the functions g{u + 1) = 
Wt^'^\wi, . . . ,Wi^-2){u + l)/ SiQ-2{u + l) and R{u) do not have common zeros. Therefore, 
the discrete Wronskian Wt^'^\wi, . . . , is divisible by R{u). □ 

3.3. Discrete exponents. For A G C, there exists an increasing sequence of nonnega- 
tive integers {ci < • ■ ■ < cm} and a basis {/i, . . . , Jm} of U such that for « = 1, . . . , M, 
we have /j(A + j) = for j = 0, . . . , q — 1 and /i(A + q) 7^ 0. This sequence of integers 
is defined uniquely and will be called the sequence of discrete exponents of U at A. We 
say that the basis {/i, . . . , /a/} agrees with exponents at A. 

We say that A is a singular point of U if the discrete exponents at A differ from the 
sequence {0 < 1 < ■ ■ ■ < M - 1}. 

For i = 1, . . . , M, introduce the local frame-type polynomials 

a-i 

Qi = n (^-^-3) 

j=Ci-i-i+2 

where cq = —1. Notice that 

• roots of each Qi are simple, 

• deg Qi = Ci- Cj_i - 1, 

• sets of roots of different polynomials do not intersect, 

• the union of roots of all Qi is the sequence A, A + 1, . . . , A + cm — M. 

Theorem 3.3. Let ci < ■ ■ ■ < cm be the discrete exponents of U at A. Then 

(i) The discrete Wronskian of U is divisible by 

M Cs-M M M-i+l 

n n = n n 

s=l j=s-M i=l j=l 

In particular, the total degree of this divisor is YlfLii^i — ^ + !)■ 

(ii) If Si, . . . , Sm o,i"g the polynomials from Lemma W^ then for any k, the polynomial 
Sk is divisible by 

k k-i+l 

n n Q^i^+J-^) ■ 

1=0 j=l 
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□ 

Corollary 3.4. Assume that the exponents of U at X have the form 

0<...<M-L-l<M-L + h<M-L + l2<---<M-L + lL 

for a suitable L, 1 < L < M , and < li < ■ ■ ■ < II. Then the local frame-type 
polynomials have the form Qi = 1 for i = 1, . . . , M — L, 

Qm-l+1 = - A - j) , QM-L+i= Yl i^-^-j) 

j=0 i=ii_i-i+2 

for i = 2, . . . , L, and the discrete Wronskian of U is divisible by 

n n ("-^-j). 

k=l j=k-L 

In particular, the total degree of this divisor is X]fc=i(^fc — ^ + !)• 

3.4. Proof of Theorem 13. 8L We shall prove part(i). Part (ii) is proved analogously. 
Let be a basis in U that agrees with exponents at A. Consider the 

matrix- valued function F{u) = [fj{u + k — l)]j.k=i,...,M ■ 

Lemma 3.5. For t E C, if the corank of F{t) is r, then the discrete Wronskian 
Wr('^)(/i, . . . , /m) is divisible by {u - ty . □ 

For i G Z, set dj = #{s < M \ j + M < cj. For j > 0, set rj = dj, and for j < 0, 
set Vj = max (0, dj + j). 

It is easy to see that for j < 0, the number rj can be also defined as < M | s < 
j + M< cj. 

Lemma 3.6. For j G Z, the corank of F{X + j) is not less than rj. 

Proof. If j > 0, then F(A+j) has dj zero rows produced by the functions fM-dj+i, ■ ■ , fn- 
Hence, the corank of F{\ + j) is at least dj = vj. 

If —M < j < 0, then the rows produced by the functions fM-dj+i, • • • , /m have zeros 
everywhere except in the first —j columns. Hence, the corank of F{X + j) is at least 
d,-{-j). □ 

By Lemma Eini the discrete Wronskian of U is divisible by Y['j=i-M{'^ " A — j)'^ • That 
can be written as n*Ii YYj''=s-Mi'^ " A — j). Theorem 13.31 is proved. 

3.5. Numerically non-degenerate space of quasi-exponentials. Let U he a space 
of quasi-exponentials like in Section 13.11 

Let Ai, . . . , A„ G C be such that Aj — Aj ^ Z for i ^ j. For i = 1, . . . , n, let the 
exponents of U at Aj have the form 

{0 < ■■ ■ < M - Ni - 1 < M - N, + na < M - Ni + n,2 < ■ ■■ < M - Ni + ni^vj 
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for a suitable iVj, 1 < A^j < M, and < riii < ■ ■ ■ < niNi- Set 

N = + ... + Nn. 
We say that the space U is a space of quasi- exponentials with data 

where a = 1, . . . , m, 6 = 1, . . . , M^, z = 1, . . . , n, j = 1, . . . , A^j. 

We say that the space U is a numerically non- degenerate space of quasi- exponentials 
with respect to data D if 

m Ma n Ni 

(3-1) 5^5^ (m,, + l-6) = + 

a=l 6=1 i=l j=l 

For 2 = 1, . . . , n, let Qa, . . . , QjAf be the local frame- type polynomials associated with 
the point Aj. For k = 1, . . . , M, define 

n 

(3.2) Qkiu) = II Q,k{u) . 

i=l 

Lemma 3.7. If U is a numerically non- degenerate space of quasi- exponentials, then the 
discrete Wronskian of U is given by the formula 

m M M-k+l 

wrj?) = n^f"" n n 

a=l k=l j=l 

Moreover, if Si, ... , Sm are the polynomials from Lemma U^ then for any i, the poly- 
nomial Si is divisible by 

i i—k+1 

n n Qki^+j-^) ■ 

k=0 j=l 

The lemma follows from Lemmas 13. 13 .21 and Theorem 13.31 

3.6. Fundamental difference operator. The monic fundamental difference operator 
of a space of quasi- exponentials U is the unique monic linear difference operator 

Du = T^' + BiT^'-' + --- + BM-iru + BM 
of order M annihilating U . Here 

• Wrg^, 
Bi = (-1)'— , 
Wr^' 

where Wr[^'* is the discrete Wronskian of a basis {/i,...,/m} of U, Wr[^''j is the 
determinant of the M x M-matrix whose j-th row is 

fjiu + 1), f^{u + N~t-l), f^(u + {N-t + l), f,iu + N) . 
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Clearly, Bi, . . . , Bm are rational functions. 

Lemma 3.8. For any i, the function Bi has a limit as u tends to infinity. Denoted this 
limit by -Bj(oo). Then 

M 

X^' + 5i(00)x'''-l + ... + 5m(00) = W{x-Zaf''^ . □ 

a=l 

Lemma 3.9. For any i = 1, . . . , M, the junction 

n Ni 

Hu) YlYliu - Xi - Uij + Ni) 

i=i j=i 

is a polynomial. 

Proof. Let A be one of the points of the set {Ai, . . . , A„}. For such a A, in the proof of 
Theorem 13.31 we defined the numbers dj and Vj for j G Z. 

For j G Z, define the new numbers Pj as follows. Set pj = dj^i for j > 0, and set 
Pj = max(0, j + dj+i) for j < 0. 

By the construction, we have pj < Vj. 

The reasons in the proof of Theorem 13.31 show that the discrete Wronskian Wi^^ is 
divisible by 

CM-M 

x{u)= n («-A-jp . 

j=l-M 

Similar reasons show that for any A;, the determinant Wr[^|, is divisible by 

CM-M-l 

Y{u)= n {u-X-jY^. 

j=l-M 

As in the end of the proof of Theorem 13. 3t we have 

M Cs-M M Cs-M 

x{u)=i[ n («-A-j)= n n (^-^-^o > 

s=lj=s-M s=l j=s-M 

Cs>S 

where in the second expression we excluded the empty products over j. Similarly, 

M Cs-l-M M Cs-l-M 

Y{u)= n n (^-^-^■)=n n (^-^-^o, 

s=l j=s-M s=l j=s-M 

where the first expression is analogous to the second expression for X[u), and the second 
expression may contain certain empty products over j. 
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Using the second expressions for X{u) and Y{u), we get 

M 

(3.3) X{u) = Y{u) Y[{u-X-Cs + M) . 

s=l 

Cs>S 

Now if A = Aj, we shall provide X{u) and Y{u) with index i. Calculating the product 
in ()3.3|) for A = Aj, we get 

Xi{u) = Yi{u) JJ (m - A, - Uij + Ni) . 

i=i 

Multiplying this formula over i = 1, . . . , n, we conclude that for any k, the product 

n Ni 

i=i j=i 

is divisible by the discrete Wronskian WT^\ □ 

Define the regularized fundamental difference operator of the space U as the linear 
difference operator 

Du = Bo{u)t^ + Bi{u)r^'~' + ... + Bm{u) 

of order M with polynomial coefficients, which annihilates U and such that its leading 
coefficient Bq is a monic polynomial of the minimal possible degree. 
We have 

deg Bq = deg Bm > deg Bi for z = 1, . . . , M — 1 , 
by Lemma (3. 8[ 

We say that the space U with data { m, Mq, niab, Za, n, Ni, riij, Aj } is a non-degenerate 
space of quasi- exponentials if U is numerically non-degenerate and 

deg 5o = N = iVi + ■ ■ ■ + iV„ . 

By LemmaEini if the space U with data { m, Ma, mat, Za, n, N^, Uij, Aj } is non-degenerate, 
then 

n Ni 

Bo = n n " ~ + • 

i=i j=i 

Example. Let U be the vector space spanned by the quasi-exponentials u and u{u — 1). 
Then M = 2, m = 1, Zi = 1, mu = 1, = 2. Let Ai = 0. The exponents of U at Ai 
are 1 and 2. Then nu = 1, = 2, N = Ni = 2. Equality p.ip takes the form: 2=2. 
We have 

Du = u{u + 1) - 2u{u + 2) r„ + (m + l){u + 2) . 
Hence U with this data is non-degenerate. 
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Example. Let U be the vector space spanned by the quasi-exponentials u and (—!)"«. 
Then M = 2, m = 2, Zi = 0, Z2 = -1, Mi = 1, rriu = 1, M2 = 1, msi = 1. Let Ai = 0. 
The exponents of U at Ai are 1 and 2. Then rin = 1, ni2 = 2, N = Ni = 2. Equahty 
()3.1|) takes the form: 2=2. Hence U with this data is numerically non-degenerate. We 
have 

Du = u tI + (m + 2) . 
Hence U with this data is degenerate. 

Example. Let U be the vector space spanned by the quasi-exponentials u and (— 
Then M = 2, m = 2, Zi = 0, Z2 = -1, Mi = 1, mn = 1, M2 = 1, msi = 1. Let 
Ai = —1. The exponents of U at Ai are 0,3. Then A'^i = 1, nu = 2, N = 1 and 
Equality p.ip is 2 = 2. Therefore, U with this data is numerically non-degenerate. 
With this data, we have 

Du = utI + (m + 2) . 
Hence, U with this data is non-degenerate. 

Theorem 3.10. Assume that the space U with data { m, M^, m^b, Za, n, Ni, riij, Xi} is a 
non- degenerate space of quasi-exponentials. Then 

(i) We have 

m n Ni 

Bm = n(-^-)'''"nn i^-^^+j) ■ 

a=l i=l j=l 

(ii) Write 

Du = n^Ao(r„)+^i^-%(r„) + --- + A;v(r„) 
where Aj(r^j) is a polynomial in with constant coefficients. Then 

(iii) The polynomials Aq, . . . , Am have no common factors of positive degree. 

Corollary 3.11. If the space U is non-degenerate with respect to a data {m, Ma,mab, 
Za, n, Ni,nij, Xi} , then the data is determined uniquely. 

The corollary follows from part (i) of the theorem. 

Proof of Theorem VJ.llA Part (ii) follows from Lemma [3. 81 Part (iii) follows from the fact 
that U does not contain exponential functions z". 
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Let Qi, . . . , be the polynomials introduced in ()3.2|1 . To prove part (i) it is enough 
to notice that 

Bm{u) _ , ^.^. Wr^'^^ln + l) _-,^ Qfc(n + M + l-fc) fr^^^ 
^ ' W.(d)/„A 11 nAn,\ 11 



fc=l ^"-^ ' a=l 

n - A, + J 



i=\ j=l a=l 



□ 



3.7. Regularized Conjugate Space. Let V with data { m, M^, mat, ^a, n, A^j, n^j, Aj } 

be a non-degenerate space of quasi-exponentials as in Section Let Wr^^"* be the dis- 
crete Wronskian of U and let Bm{u) be the last coefficient of the regularized fundamental 
difference operator of U. 

The complex vector space spanned by all functions of the form 

with fi E V has dimension M. This space is denoted by and called regularized 
conjugate to U . 

Lemma 3.12. For any g G f/"'", the function 

i=l j=0 

is holomorphic in C. 

Proof. Let Qi, ■ ■ ■ ,Qm be the polynomials introduced in (j3.2j) . Let g he a. function in 
(HI. Then r„(Wr('i)(/i, . . . , /m-i)) is divisible by 

A/-1 A/-2 

n + ■ n Q^^^ +')■■■ qm-i{u + 1) 

i=l i=l 

while 

m M M-1 

a=l j=l j=l 

Hence, the possible poles of g come from the product Bm{u)Qi{u) . . .Qm{u) which 
remains in the denominator of g. But this product is exactly the product in Lemma 
EH □ 

For every i = 1, . . . , n and j = 1, . . . ,Ni, fix a function gij in U which is equal to 
zero at M = A^, Aj + 1, . . . , Ai -|- M — iVj + n^j — 1 and which is not equal to zero at 
u = Xi + M - Ni + riij. 
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Lemma 3.13. For given i = 1, . . . ,n, j = 0, . . . , Ni — 1, let fi, . . . , /m-i be a collection 
of functions in U containing the functions gij+i, ■ ■ ■ i9i,Ni cind let 

p,. ^ ru{WT^'\fl,---,fM-l)) 

Then for j = 0, the function F has no poles at 

u = Xi- N„ Xi- Ni + 1, . . ., Xi- Ni + UiNi . 

For j > 0, the function F has no poles at 

u = Xi- Ni + Uij + 1, Xi- Ni + Uij + 2, . . . , Xi- Ni + Ui^^ . 

Moreover, if fi, . . . , /a/_i is a generic collection of functions in U containing the func- 
tions gi,j+i, gi,j+2, ■ ■ ■ , 9i,Ni, then the function F has a nonzero residue atu = Xi-Ni+Uij. 

Proof. The first two statements of tlie lemma are proved in tlie same way as Lemma 

We sliall prove tliat tlie residue of F at u = Xi — Ni + riij is nonzero first assuming 
that Uij > Ni. To prove that the residue is nonzero it is enough to show that 

(3.5) ord„=A,-7v,+n,, Wr[?^ = 1 + ord„=A,-7v,+n,,+iWr[?^ = Ni - j + 1 , 
and 

(3.6) ord„=A,-jv,+n.,+iWr[?^ = ord„=A,-jv,+n,,+iWr('^) (/i, • • • , /m-i) • 

Here ord„=A/ denotes the order of zero of the function / at u = A. 

Equalities (j3.5|) follow from the numerical non-degeneracy of the space U and the 
condition Uij > Ni. 

Since the collection fi, . . . , fu-i is a generic collection containing the functions gi,j+i, 
. . . , gi^Ni, the discrete Wronskian Wr'''^-*(/i, . . . , fM-i,gij) is nonzero and proportional to 
Wrl^"*. Expanding the determinant Wr'''^-*(/i, . . . , fM-i,gij) with respect to the last row, 
we have 

(3.7) Wt'^\u) = const {g.i,{u + M - 1) Wi^'^^fi, . . . , /a/-i)(m) - . . . 

- (-l)^^^7,,(..)r„(Wr(^)(A,...,/M_i))(«)) . 

where const ^ 0. 

The order of Wr*-'^-'(/i, . . . , fM-i) at u = Xi — Ni + nij + 1 is not less than Ni — j. This 
follows from Theorem 13 . 31 applied to the functions /i, . . . , fu-i- A similar reason shows 
that the order at m = A,; — Ni+ Uij + 1 of all of the other (M — 1) x (Af — 1) minors in 
the right hand side of ()3.7|) is also not less than Ni — j. 

By the construction, the function gij is nonzero at u = Xi — Ni + Uij and is zero 
at M = Aj — A'j + riij — I for / = 1,...,M — 1. Therefore, the only term in the right 
hand side of equahty ()3.7|) that can have order Ni — j at u = Xi — Ni + riij + 1 is 
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9ij[U 



+ M-1) Wr('^)(/i, . . . , /m-i)(m). Since ord,= 



■Xi-Ni+nij + 



iWr 



(d) 
U 



Ni— j, this shows 



that the orders of Wrl^-* and Wr*-'^''(/i, . . . , /a/-i) at u = Xi — + equaL 

To prove that the residue oi F aX u = Xi — Ni + rtij is nonzero in the case riij < Ni/ii 



is enough to show that 



ord„=A,-iv,+n„Wr 



(d) 

u 



ord„=A,-7v,+n„+iWr 



(d) 

u 



Uij -j + l 



and 

ord„=A,-7v,+n,,,+iWr[?^ = ord„=A,-7v,+n,,+iWr('^)(/i, . . . , /m-i) 
The proof is similar to the proof of equahties ()3.5|) and ()3.6|) . 



□ 



Theorem 3.14. // Du = Bq{u)t^^ + Bi{u)t^^-^ + • ■ ■ + TuBm-i{u) + Bm{u) is the 
regularized fundamental difference operator of U, then the operator 

D\j = t!:'Bm{u)+t^'-^Bm-i{u) + --- + tIB^{u) + B^{u) 

annihilates . 

Proof. Consider the scalar equation -^^Duy = with respect to an unknown function 

y{u). For i = 1, . . . , M, introduce Wi = T^~^y and present the equation as a system of 
first order equations 



TuWm 



Biju) 
' B.iu 



-wm - 



Bm{u) 
B.(u) 



Wi 



TuWi = Wi+i 



for 2 = 1, . . . , M — 1. For the column vector w = {wi, . . . , wm), the system can be 
presented as a matrix equation TuW = Cw with the M x M-matrix 



C 








Bm{u) 
Bo{u) 



1 




-Ba/-i(m) 
Bo{u) 





1 



Bm-2(u) 
Boiu) 








B2{u) 
'Boiu) 








1 

Bo{u) / 



Let \& be a fundamental M x M-matrix of solutions, = C\E'. Then r„\E' 
where 



-1 



I Bm-i(u) 
Bm{u) 
1 





Bm-2{u) 

Bm{u) 



1 



Bm~3(u) 

Bm{u) 






Biju) 
Bm{u) 






Boju) 
Bm(u) 
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For a row vector (t>i, . . . , vm) the equation TuV = vC~^ has the form: 

B B 

TuVm = -—vi, TuVm-1 = --^vi + Vm, 

ruVM-2 = -^—Vl + Vm-1, • • • , TuVi = Vi + V2 

DM -DM 

This system reduces to the scalar equation D^jVi = where 



Bm-i{u) , , fii(n) , Bo{u) 
^ Bm{u) ^'"^ ""BMiu)^ Bm{u) • 

Thus the kernel of the difference operator consists of the first row entries of the 
matrix \E'~^. 

If {/i, . . . , /m} is a basis of U, then 

/ /iH /2(n) fM{u) \ 

^_ fliu + l) f2iu+l) fM{u + l) 



\Mu + M-l) f,{u + M-l) f,,^(u + M-l) J 

The formula for the inverse matrix elements shows that DIj annihilates the functions of 
the form r„(Wr^'^-*(/i, . . . , fM-l))/^'c^\u). Then the operator Dfj annihilates the func- 
tions of the form ()3.4|) . since Dfj = ■ Bm where -Bm is the operator of multiplication 
by the function Bm- □ 

4. Integral Transforms 

4.1. Mellin-type transform. Let be a non-degenerate space of quasi-polynomials 
with data S}y = {n, Ni, riij, Aj, m, Ma, rriab, Za} where i = 1, . . . ,n, j = 1, . . . , Ni, a = 
1, . . . , m, 6 = 1, . . . , Ma- Let V'^ be the space regularized conjugate to V. 

For a = 1, . . . , m, denote by 7^ a small circle around Za in C oriented counterclockwise. 

Denote by U the complex vector space spanned by functions of the form 



(4.1) Uu) = / x-fix)dx, 

J-ya 

where a = 1, . . . , m, f eV'^ . The vector space U is called bispectral dual to V. 
Theorem 4.1. Let V be a non- degenerate space of quasi-polynomials with data 

Dv = {n,Ni,nij,\i, m,Ma,mab,Za} 
where i = 1, . . . ,n, j = 1, . . . , Ni, a = 1, . . . ,m, b = 1, . . . , Ma- Then 

(i) The space U is a non- degenerate space of quasi- exponentials with data 

= {Tn,Ma,Tnab,Za, n,Ni,nij,\i + Ni} 
where a = 1, . . . ,m, b = 1, . . . , Ma, i = I, . . . ,n, j = 1, . . . , Ni. 
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(ii) Let Dy = X]f=i XljLi regularized fundamental differential op- 

erator of V where Aij are suitable complex numbers. Then 

M N 

i=i j=i 

is the regularized fundamental difference operator ofU. 
The theorem is proved in Section 

4.2. Fourier-type transform. Let U he a non-degenerate space of quasi-exponentials 
with data Du = {m, Ma, mat, Za, n, N^, riij, Aj} where a = 1, . . . ,m, b = 1, . . . , Ma, i = 
1, . . . ,n, j = 1, . . . , Ni. Let be the space regularized conjugate to U. 
For i = 1, . . . ,n, consider the arithmetic sequence 

6, = {K~N„Xi-Ni + l,Xi-N, + 2,...,Xi-Ni + niN,} 

consisting of ni^. + 1 terms. 

For i = 1, . . . ,n, fix a non-self-intersecting closed connected curve 7^ in C oriented 
counterclockwise and such that the points of the sequence &i are inside 7^ and the points 
of other sequences &j for j i are outside 7^. 

Denote by V the complex vector space spanned by functions of the form 

(4.2) fi{x) = [ x-f{u)du, 

where i = 1, . . . ,n, f eU^ . The vector space V is called bispectral dual to U. 
Theorem 4.2. Let U be a non-degenerate space of quasi-exponentials with data 

Du = {m,Ma,mab,Za, n,Ni,nij,Xi} 

where a = 1, . . . ,m, b = 1, . . . , Ma, i = 1, . . . ,n, j = 1, . . . , Ni. Then 

(i) The space V is a non-degenerate space of quasi-polynomials with data 

Dv = {n,Ni,nij,Xi- Ni, m,Ma,mab,Za} 

where i = 1, . . . ,n, j = 1, . . . , Ni, a = 1, . . . ,m, b = 1, . . . , Ma. 

(ii) Let Djj = X^ili Sjli ^ij'^^^u ^6 ^he regularized fundamental difference operator 
of U where Aij are suitable complex numbers. Then 

N M 
i=l j=l 

is the regularized fundamental differential operator ofV. 
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The theorem is proved in Section 031 

Theorems 14.11 and 14.21 imply that if is a non-degenerate space of quasi-polynomials 
and U is bispectral dual to V, then V is bispectral dual to U. Similarly, if ?7 is a 
non-degenerate space of quasi-exponentials and V is bispectral dual to U, then U is 
bispectral dual to U. 

4.3. Proof of Theorem 14. IL The exponents of V"^ at Za are 

{-ma,Ma - 1< • ■ ■ < -rriai -l<0<l<---<N-Ma-l}. 

Integral ()5.H) is nonzero only if / has a pole aX x = Za- If / has a pole at x = Za of 
order —rriab — 1, then integral ()5.1|) has the form z^qab{u) where qab is a polynomial in 
u of degree rriab- Thus f/ is a space of quasi-exponentials of dimension M generated by 
quasi-exponentials z^qab{u) where a = 1, . . . ,m, 6=1,..., mab and qab is a polynomial 
in u of degree niab- 

It is clear that the operator D"^ = Xlj^i ^ij ^"'^u annihilates U. Write 

= Bq{u)t^' + ■ ■ ■ + Bm-i{u)tu + Bm{u) 

where Ba{u) are polynomials in u with constant coefficients. Lemma 12.31 implies that 
Bq{u) = 11^=1 ~ -^j)^^ the polynomials Bq, . . . ,Bm have no common factor of 
positive degree. Therefore, is the regularized fundamental operator of U. 

The functions x^^pij{x), i = 1, . . . , n, j = 1, . . . ,Ni, form a basis of V. For given 
order all basis functions of V except the function x^^pij{x). Denote by Wrjj the 
Wronskian of this ordered set of iV — 1 functions. The functions 

f ^ "^Tjjjx) 

^''^ ' mv{x)x^ YlZA^-^a)''^ 

form a basis in V'^ . Such a function fij has the form x~^^rij where rational 
function in x. We have 

ord^=o Tij = -Ni , orda;=oo nj = -Uij - M -1 . 

Consider the following element of U, 

m „ 

a=l >^U-i7a 

If M = Aj + m where Ni < m < Uij + M — 1, then the integrand is a rational function with 
zero residues at and oo. Hence Fij is zero at m = Aj + m for m from this arithmetic 
sequence. This remark together with Theorem 13.31 proves that U is a non-degenerate 
space of quasi-exponentials with data { m. Ma, niab, Za, n, Ni, Uij, Aj + Ni}. □ 
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4.4. Proof of Theorem 14. 2L By Lemmas 13.121 and for any i = 1, . . . ,n, the 

functions fi{x) in ()4.2|) have the form x'^^~^^Pij{x) where j = l,...,iVj and Pij is a 
polynomial of degree not greater than riij. Moreover, there exists a function x^^~^^pij{x) 
with Pij of degree exactly equal to riij. 

The functions z^qab{u), a = 1, . . . ,m, 6=1,..., Ma, form a basis of [/. For given a, 6, 
order all basis functions of U except the function z^qabiu) and denote them /i, . . . , /m-i- 
The corresponding function 

, . . r„(WrW(A,...,/M-i)) ^ 
5M(u)Wr[?)(n) 

has the form z~'^rab{u) where Tab is a rational function in u and ord„=oo ^afe = Ma — 
^ab — N — 1. Consider the following element of V, 



n 



Fab{x) = ^ifab)i{x) = / X'^fij{u)du. 
i=l 

li X = Za, then the integrand is a rational function in u which tends to zero as u tends 
to infinity. Denote by the i-th derivative. Then F^l\za) = for i = 0,1, . . . , N — 

Ma + mab - 1, and F^f J ^ 0. 
This reason proves part (i) of the theorem. 

From Theorem 13.141 and formulas for the Fourier-type integral transform, it follows 
that the differential operator Xlili X]j=i x^i^xdxY annihilates V. From part (iii) of 
Theorem 13.101 it follows that this operator is the regularized fundamental differential 
operator of V. □ 



5. Special spaces 

5.1. Spaces of quasi-polynomials of an (n. A, m, z)-type. Let be a natural num- 
ber, > 1. Let n = {rii, . . . , hn) be a vector of nonnegative integers. 

Let A = (Ai, . . . , Aat) G C^. Assume that Aj — Aj ^ Z for i ^ j. 

For i = 1,. . . ,N, let Pi G C[x] be a polynomial of degree rii such that Pi{0) ^ 0. 
Denote by V the complex vector space spanned by functions x^^Piix), % = 1,. . . ,N. 
The dimension of V is A^. 

Let z = {zi, . . . , zm), M > 1, be a subset in C containing all singular points of V . 
Assume that for a = 1, . . . , M, the set of exponents of V at Za has the form 

{0<1<---<A^-2<A^-1 + ma} , ma > . 

We have 

i=l ni = La=l ■ 

We call the pair (y,z) a space of the {n, X,m, z)-type or a special space of quasi- 
polynomials. 
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Let Dy be the monic fundamental differential operator of V . The operator 



DV = Ua=l - ^a) Dv 

is called the special fundamental differential operator of the special space {V,z). 
Write the special fundamental differential operator in the form 

Dv = Ao(x)(a;9,.)^+ A,{x){xd,f-' + A2{x){xd,f-^ + ■■■ + Aj,{x). 

Then by Lemma l2.3[ all of the coefficients Aq, . . . ,A]sf are polynomials in x of degree 
not greater than M. If we write 

Dv = x^'Boixd,) + x''-^B^{xd^) + ■■■ + xBM-i{xd,) + BMixd,), 
where Bi is a polynomial in xd^ with constant coefficients, then 

TV M N 

Bo = H {xd, - A, - m) , Bm = i-l f' H^aH {xd, - A,) . 

i=l a=l i=l 

5.2. Spaces of quasi-exponentials of an (m, z, n, A)-type. Let M be a natural 
number, M > 1. Let m = (mi, . . . ,mM) be a vector of nonnegative integers. 

Let z = {zi, . . . , Zm) be distinct nozero complex numbers with fixed arguments. 

For a = 1,...,M, let qa G C[u] be a polynomial of degree m^. Denote by U the 
complex vector space spanned by the functions z'^qa{u), a = 1, . . . , M. The dimension 
of U is M. 

Assume that for some natural number > 1, there exists a subset of distinct numbers 
A = (Ai, . . . , Aat) in C with three properties: 

• Xi — \j ^ Z for i 7^ j. 

• For i = 1, . . . , N, the set of discrete exponents of U at Aj has the form 

{0<1 < ■■■ < M -2< M -l + m} , ni>0 . 

EN v-^Af 
i=l ni = La=l 

In this case we call the pair {U, A) a space of quasi-exponentials of the (m, z, n, \)-type 
or a special space of quasi-exponentials. 

Let Djj be the monic fundamental difference operator of U . The operator 

Du = nil iu-\i-ni + 1) Dv 

is called the special fundamental difference operator of the special space {U,X). 
Write the special fundamental difference operator in the form 

Du = Bo{u)t^'+ B,{uyj''+B,{uyj-^+ ■ ■ ■ + Bm{u). 

Then by Lemma and Theorem l3.10l all of the coefficients i?o, • • • , Bm are polynomials 
in u of degree not greater than and 

M N 

Bm = {-lf'\[za\[{u-K + l). 

a=l i=l 
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If we write 

Djj = U^Ao{Tu) + U^-^Ai{Tu) + ■■■ + uAN-liTu) + ANiTu), 

where v4j(r„) is a polynomial in r„ with constant coefficients, then 

M^u) = Hall - ^a) ■ 

5.3. Special Mellin-type transform. Let {V, z) be a space of quasi-polynomials of 
an (n, A, m, 2)-type. Let V* be the space conjugate to V and defined in Section YIM 

For a = 1, . . . , M, let 7^ be a small circle around Za in C oriented counterclockwise. 
Denote by U the complex vector space spanned by the functions of the form 

„ M 
(5.1) fa{u) = / X''f{x)x-''\[{X-Zar^dx, 

"^T" a=l 

where a = 1, . . . , M, / G 

Theorem 5.1. Lei (V, 2) 6e a sj»ace of quasi-polynomials of an 
(n, (Ai, . . . , Aat), m, (zi, . . . , zm))- i?/l>e. Then 

(i) r/ie sj>ace U is a space of quasi- exponentials of the 
(m, (^1, . . . , ztv), n, (Ai + 1, . . . , Atv + l))-type. 

(ii) Lei = Xlfii XljLi ^ij^^^dxY be the special regularized fundamental differen- 
tial operator of V where Aij are suitable complex numbers. Then 

M N 

E E 

i=i j=i 

is the special regularized fundamental difference operator of U. 
The proof is similar to the proof of Theorem 14.11 

The special space of quasi-exponentials {U, A) is called special bispectral dual to the 
special space of quasi-polynomials (y,z). 

5.4. Special Fourier-type transform. Let {U, A) be a space of quasi-exponentials of 
an (m, 2, n, A)-type. Let Wr^^-* be the discrete Wronskian of U and let Bm{u) be the 
last coefficient of the special regularized fundamental difference operator of U. 

The complex vector space spanned by all functions of the form 

.52) r.(WrW(/i,...,/M-i)) 

with fi & V has dimension M. This space is denoted by U' and called special regularized 
conjugate to U . 

For i = 1, . . . , A^, consider the arithmetic sequence 

&i = { Ai - 1, Ai, Ai + 1, . . . , Aj + rij - 1 } 

consisting of + 1 terms. 
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For i = 1, . . . , A^, fix a non-selfintersecting closed connected curve 7^ in C oriented 
counterclockwise such that it encircles the sequence &i and does not contain inside or 
intersect with points of other sequences &j for j 7^ i. 

Denote by V the complex vector space spanned by functions of the form 



where z = 1, . . . , n, f G U' . The vector space V is called special bispectral dual to U. 

Theorem 5.2. Let {U,X) be a space of quasi- exponentials of an 
{m, (zi, ...,ZM),n, (Ai, . . . , \N))-type. Then 

(i) The space U is a space of quasi- exponentials of the 
(n, (Ai - 1, . . . , Aat - 1), m, (zi, . . . , ZN))-type. 

(ii) Let Du = XljLi ^ij'^^'^l special regularized fundamental difference 
operator of U where Aij are suitable complex numbers. Then 



is the special regularized fundamental differential operator of V. 
The proof is similar to the proof of Theorem 14.21 

The special space of quasi-polynomials (V, z) is called special bispectral dual to the 
special space of quasi-exponentials (f/, A). 

Theorems 15.11 and 15.21 imply that if is a special space of quasi-polynomials and U 
is special bispectral dual to V, then V is special bispectral dual to U. Similarly, if U 
is a special space of quasi-exponentials and V is special bispectral dual to U, then U is 
bispectral dual to U. 



6. Special spaces and solutions of the Bethe ansatz equations 

6.1. Critical points of master functions and special spaces of quasi-polynomials. 

Let (y, z) be a space of quasi-polynomials of an (n, A, m, 2;)-type. We construct the as- 
sociated master function as follows. Set 




M N 



i=i 3=1 



Hi = Ui+i H h riAT , 



z = 



l,...,iV-l . 



Consider the new ni -(-■•• -(- tt-tv-i auxiliary variables 




1 ) • • • ) ''n2 1 



(2) ,(2) 



t 



(N-l) 



1 



5 • • • ) ''njv-i / • 
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Define the master function 

M N Ui 

(6.1) $(t<">;A;m;z) = JJ ;2™4Ai+m./2) J-j- (^^ _ ^,)-»-^ JJ JJ 

a=l l<a<6<M i=l j=l 

M ni Af-1 Af-2 ni "i+i 

X nn(4" nncf -4'')^nnn(*f-r')-'- 

a=l j=l 1=1 j<j' 1=1 j=lj'=l 

The master function is symmetric witli respect to tlie group S^^ = Sjj^ x ■ ■ ■ x Sfijv-i of 
permutations of variables t^*'' preserving the upper index. 

A point t^"^ with complex coordinates is called a critical point of $( ■ ; A; m; z) if the 
following system of ni + ■ ■ ■ + n^^i equations is satisfied 

X — X 1 m 2 "2 

(6-2) 7(1^ + XI ,(1) _ ^ 5Z ,(1) _ ,(1) + ,(1) _ ,(2) 







Atv-I — Atv — i ^ 



2 1_ 



,(A^-l) Z^ ,(iV-l) _ ,(7V-1) Z^ AN-l) _ ,(7V-2) 







where j = 1, . . . , fii in the first group of equations, i = 2, . . . , — 2 and j = 1, . . . , in 
the second group of equations, j = 1, . . . , n^_i in the last group of equations. 
In other words, a point i^"'^ is a critical point if 

^ = °' ^ = 1,---,A'-1, J = l,...,n.. 

In the Gaudin model, equations ()6.2|1 are called t/ie Gaudin Bethe ansatz equations. 

The Sfi-orbit of a point t^"'^ G C"'^"'" ""'""^-! is uniquely determined by the N — 1-tuple 
?/*^"* = (y^, . . . , ?/7v-i) of polynomials in x, where 

= n-=i(^-f ), ^ = i,...,iv-i. 

We say that y represents the orbit. Each polynomial of the tuple is considered up to 
multiplication by a nonzero number since we are interested in the roots of the polynomial 
only. 

We say that t^"^ G C"'i+'"+"'^-i is Gaudin admissible if the value <l>(t^"^; A; m; z) is 
well defined and is not zero. 

A point t^"^^ is Gaudin admissible if and only if the associated tuple has the following 
properties. 

• For a = 1, . . . , M, if rua > 0, then yi{za) ^ 0. 
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• For all i, yi{0) ^ 0. 

• For all i, the polynomial yi has no multiple roots and no common roots with yi_i 
or yi+i. 

Such tuples are called Gaudin admissible. 
Return to (V, 2), a space of an (n, A, m, z)-type. 

The space V = {x^^pi{x), . . . , x'^'^p^ix)) determines the A^— 1-tuple = {yi, . . . , yN-i] 
of polynomials in x, where 

yN-i =Pn, y^ = ^('^-"(f-'"''-^'+i--^ivWr . . . , x^^pn{x)) 

for 2 = 1, . . . , iV - 2. 

We call the special space (V^, z) Gaudin admissible if the tuple is Gaudin admissible. 

Theorem 6.1. |MVniMVI] 

(i) Assume that the special space (V, z) is Gaudin admissible. Then the tuple 
represents the orbit of a critical point of the master function. 

(ii) Assume that t^'^^ is Gaudin admissible and t^"^ is a critical point of the master 
function. Let y = (yi, . . . ,yN-i) be the tuple representing the orbit oft^'^K Then 
the differential operator 



yi J \ y2 

yN-i J 

of order N is the monic fundamental differential operator of a Gaudin admissible 
special space of quasi-polynomials (V, z) of the (n, A, m, z)-type. 
(iii) The correspondence between Gaudin admissible special spaces of quasi-polynomials 
of the (n, A, m, z)-type and orbits of Gaudin admissible critical points of the mas- 
ter function described in parts (i), (ii) is reflexive. 

This theorem establishes a one-to-one correspondence between Gaudin admissible spe- 
cial spaces of quasi-polynomials of the (n, A, m, z)-type and orbits of Gaudin admissible 
critical points of the master function. 

6.2. Solutions of the Bethe ansatz equations and special spaces of quasi-expo- 
nentials. Let (?7, A) be a space of quasi-exponentials of an (m, z, n, A)-type. We define 
the associated system of Bethe ansatz equations as follows. Set 

rha = rua+i H h tum , a = 1, . . . ,M - 1 . 

Consider the new rhi + ■ ■ ■ + rfiM-i auxiliary variables 

+ M _ ^.(1) ^(1) ^(2) ,(2) ,(M-1) ,(M-l)s 

L yi^ ) • • • ; t-mi) ''1 ) • • • 5 ''m2) • • • 5 ''1 ) • • • ) '^niM-i J ' 
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The Bethe ansatz equations is the following system of mi + ■ ■ ■ + friM-i equations: 

^ .(1) X ^(1) .(1) 1 '^2 .(1) .(2) , . 



fha-l (a) ,(a-l) ,(a) ,(a) . rha+i (a) (a+1) -, 

n '^b ~ '^b' TT '^b ~ '^fc' ~ TT '^6 ~ '^6' ^ Za+l 

,(a) ,{a-l) T 11 ,(a) ,(a) , . 11 



J. J. ,(a) _ _ 1 ii _ ,(a) ,1 ii _ .(a+1) 

b'=l % ^b' ^ b'=ib ^b ^b' "T b'=l ^b ^b' 



Z„ 



rhM-2 (N-1) AN-2) AN-l) (N-l) . 

n ^b ~ '^b> TT '^b ~^b' ~ ^ 

(iV-1) _ (iV-2) _ 11 (iV-1) _ (7V-1) 

fe'=l ^-b ''b' bV'' ^' b' + ^ ^ 

where 6 = 1, . . . , mi in the first group of equations, a = 2, . . . , M — 2 and 6 = 1,..., fha 
in the second group of equations, 6 = 1,... ,rhM-i in the last group of equations. 

The system of Bethe ansatz equations is symmetric with respect to the group 
Tifn = X ■ ■ ■ X SmA/_i of permutations of variables t^"'^ preserving the upper in- 
dex. 

In the XXX model, equations ()6.3|) are called the Bethe ansatz equations. 

The Syji-orbit of a point t^'"^ G C^^+'-'+^^m-i jg uniquely determined by the M - 1- 
tuple y*^"*^ = (?/i, . . . , Vm-i) of polynomials in u, where 

Va = nr=i ("-4'^^) , a = i,...,M-i. 

We say that y represents the orbit. Each polynomial of the tuple is considered up to 
multiplication by a nonzero number. 

We say that G C™i+-+'^a/-i jg XXX admissible of A-type if 

,(a) / ,{a) ,{a) , Aa) , Aa) , Aa+1) Aa) , Aa+l) _ . 

^b I ''b' ' ''b 7^ ''b' ^ 5 ''b / ''b' ' ''b T ^b' ' 

for all a, 6, 6', i, and r = 0, . . . , nj. 

If t^"^^ G C™i+ "+'"*^-i is XXX admissible, then the corresponding tuple 2/*^""^ is called 
XXX admissible. 

Return to (f/. A), a special space of quasi-exponentials of an (m, z, n, A)-type. 
The space [/ = {z^qiiu)^ . . . , z'^qMiu)) determines the M— 1-tuple = {yi, . . . , 
of polynomials in u, where 

M 

VM-i = qM , = n ^fc'" ^^^'^^ (^^+1 ^"+1 (")'•••' 4f9Af (m) ) 

b=a+l 

for a = l,...,M-2. 

We call the special space {U, A) XXX admissible if the tuple y^ is XXX admissible. 

Theorem 6.2. pVTllMVi] 
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(i) Assume that the special space {U, A) is XXX admissible. Then the tuple 
represents the orbit of a solution of the XXX Bethe ansatz equations \6. "J\] . 

(ii) Assume that t^"*^ is XXX admissible and t^"^^ is a solution of the XXX Bethe 
ansatz equations \6.'^ Let y = {yi, . . . be the tuple representing the orbit 
of t^'^\ Then the difference operator 

n ( TT u-\i + l \[ yi{u + l) y2{u) 



yi{u + l) fj^u- Xi-rii + l J\ yi{u) y2{u + 1 



Z2 



?/Af_2(M + l) yM~i{u) \( yM-i{u+l) 

Tu -T-. 7 — —rr zm-1 ) \Tu J—. — Zm 

yM~2{u) yM-i{u + l) J\ yM-i[u) 

of order M is the monic fundamental difference operator of a special space of 
quasi- exponentials {U, A) of the (m, z, n, X)-type. 
(iii) The correspondence between XXX admissible special spaces of quasi- exponentials 
of the (m, z, n, \)-type and orbits of XXX admissible solutions of the XXX Bethe 
ansatz equations described in parts (i), (ii) is reflexive. 

This theorem estabhshes a one-to-one correspondence between XXX admissible spe- 
cial spaces of quasi-exponentials of the (m, z, n, A)-type and orbits of XXX admissible 
solutions of the XXX Bethe ansatz equations. 

7. FiNITENESS OF SOLUTIONS OF THE BETHE ANSATZ EQUATIONS 

7.1. Finiteness of the Gaudin admissible critical points. 

Lemma 7.1. For given fixed m, z and generic X, the master function A; m; z), 

defined in h6. has only finitely many Gaudin admissible critical points. 



The lemma follows from Lemma 2.1 in |MVlj . 



7.2. The number of orbits of the Gaudin admissible critical points. Let Ai, . . . , 

Aat G C be distinct numbers such that Aj — Aj ^ Z for i ^ j. Set n = ni -\- ■ ■ ■ -\- n^. 
Consider the complex vector space X spanned by functions x^^~^\ i = 1, . . . , N, j = 
0, . . . , rij. The space X is of dimension n -\- N. 
For z G C*, define a complete flag F{z) in X, 

F{z) = {0 = Fo{z)cF,{z)c---CEr,+N{z)=X} , 

where Fk{z) consists of all / G X which have zero at z of order not less that n + N — k. 
The subspace F^^z) has dimension k. 

Define two complete flags of X at infinity. 

Say that x^^^^ <i x'^^'^^ ii i < i' ot i = i' and j < j'. Set 

F(ooi) = { = Fo(cx)i) C Fi(ooi) C ■ ■ ■ C Fs+jv(ooi) = X } , 
where Ffc(ooi) is spanned by k smallest elements with respect to <i. 
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Say that x'^^'^^ <2 x'^*''^^' ii i > i' or i = i' and j < j'. Set 

F(oo2) = {0 = Fo(oo2) C Fi(oo2) C • • • C F^,+jv(oo2) = X} , 
where ^^(002) is spanned by k smallest elements with respect to <2- 

Denote by Gr{X, N) the Grassmannian manifold of iV-dimensional vector subspaces 
of X. Let F be a complete flag of X, 

F = {0 = Fo C Fi C ... C Fn+N = X} . 

A ramification sequence is a sequence (ci, . . . , Cat) G such that n > Ci > • • • > Cjv > 0. 
For a ramification sequence c = (ci, . . . , c^) define the Schubert cell 

ni{F) = {V e Gr{X, N) I dim{V f] F„) = i, 

n + i-ce < u < n + £+l- q+i, i = Q,...,N} , 

where cq = n, c^+i = 0. The cell Q^{F) is a smooth connected variety. The closure of 
is denoted by Qc{F). The codimension of ^ci^) 

\c\ ^ Ci + C2 + . . . + Cn . 

Every A'"-dimensional vector subspace of X belongs to a unique Schubert cell 

For a = 1, . . . , M, define the ramification sequence 

c(a) = {rua, 0, . . . , 0) . 

Define the ramification sequences 

c(ooi) = {n2-\ \-nN,n3-\ \- un, ■ ■ ■ ,nN,0) , 

0(002) = (m H h nN-i,ni H h nN-2, . . . , 0) . 

Lemma 7.2. 

• We have 

M 

J2 codimf}^(„)(F(zJ) + codimQ°(^^)(F(ooi)) + 

a=l 

codim f^^(^^)(F(oo2)) = dim Gr{X,N) = Nn . 

• Let V G Gr(X, iV). The pair {V, z) is a space of the {n, X,m, z)-type, if and 
only if V belongs to the intersection of M + 2 Schubert cells 

Q^,(i)(F(zO)nQ°(2)(F(z2))n... 

n q°(^)(f(zm)) n Q°(^^)(F(ooi)) n ni^^^^{F{oo2)) ■ 

□ 
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According to Schubert calculus, the multiplicity of the intersection of Schubert cycles 

fie(l)(i^(^l))nfie(2)(i^(^2))n... 

(7.1) n n^(^M){FizM)) n n^^^,){F{ooi)) n n,^^,){F{oo2)) 

can be expressed in representation-theoretic terms as follows. 

For a ramification sequence c denote by lI^^ the finite dimensional irreducible 
0[^-module with highest weight c. Any gl^-module L^^^ has a natural structure of 
an slTv-module denoted by L^'^K By jFu], the multiplicity of the intersection in (|7.ip is 
equal to the multiplicity of the trivial s[7v-niodule in the tensor product of slAr-modules 

(7.2) m ® . . . ® ® z<7> , ® , . 

\ ' c(l) c(M) c(ooi) c(oo2) 

Proposition 7.3 ( jMTV2j ) . The multiplicity of the trivial 51n -Tnodule in the tensor 
product (|y.ij| ) is equal to the dimension of the weight subspace of weight [ni, . . . ,71^] in 
the tensor product of qI^ -modules 

(7.3) Li'^l ® . . . Li%^ . 

Corollary 7.4. For generic X, the number of orbits of the Gaudin admissible critical 
points of the master function $(t^"'^; A; m; z), is not greater than the dimension of the 

weight space {L^c{i) ® ■ ■ ■ ® -^i^/)) ["-i' • • • ' ^n] ■ □ 

7.3. Finiteness of solutions of Bethe ansatz equations (j6.3|) . Let z = (zi, . . . , zm) 

G C^^. Let (mi, . . . ,mM-i) be a collection of nonnegative integers. We say that z is 



separating with respect to (mi, . . . , niM-i) if 

M-l , . 

n (— ) ^1 



a=l ^ ' 

for all sets of integers {ci, . . . , cm-i} such that < Ca < ma, > 0. 

Clearly, for given (mi, . . . , mM-i), a generic z is separating. 

Lemma 7.5. Let n, X, {rfii, . . . ,ThM-i) be fixed. Let z be separating with respect to 
(mi, . . . ,mj\/_i). Then the Bethe ansatz equations \6. 'J\) have only finitely many XXX 
admissible solutions of X-type. 

Proof. If the algebraic set of XXX admissible solutions of ()(j.3j) is infinite, then it is 
unbounded. Suppose that we have a sequence of solutions which is unbounded. Without 
loss of generahty, we assume that t^^^ tends to infinity for a = 1, . . . , M — 1, 6 = 1, . . . , c^, 
and remains bounded for all other values of a, h. Multiply all of the equations in ()6.3|) 
corresponding to a variables t'f^ with a = 1, . . . , M — 1, 6 = 1, . . . , c^, and pass in the 
product to the limit along our sequence of solutions. Then the resulting equation is 
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This equation contradicts to our assumption. □ 

7.4. The number of orbits of the XXX admissible solutions. Let Zi, . . . , zm G C 

be distinct numbers with fixed argument. Set m = mi + - ■ ■ + mM- Consider the complex 
vector space Y spanned by functions z^u'^, a = 1, . . . , M, 6 = 0,..., rria- The space Y 
is of dimension fh + M. 

For A e C, define a complete fiag -F(A) in y, 

F(A) = {0 = Fo(A)cFi(A)c---cF^+M(A)=y}, 

where -Fjt(A) consists of all / e y which are divisible by HJlJi^ '' {u — X — j + 1). The 
subspace Ffc(A) has dimension k. 

Define two complete flags of Y at infinity. 

Say that z'^u'^ <i z'^,u'^ if a < a' or a = a' and b < b'. Set 

F(ooi) = { = Fo(ooi) C Fi(ooi) C • • • C F^+m(ooi) = F } , 

where Ffc(ooi) is spanned by k smallest elements with respect to <i. 
Say that z'^u'^ <2 z'^,u^' if a > a' or a — a' and b < b'. Set 

F(oo2) = {0 = Fo(oo2) C Fi(oo2) C • • • C F^+mM ^Y} , 

where ^^(002) is spanned by k smallest elements with respect to <2. 

Denote by Gr{Y, M) the Grassmannian manifold of M-dimensional vector subspaces 
of Y. Let be a complete flag of Y, 

F = { = Fo C Fi C ... C F„-,+M = Y} . 

For a ramiflcation sequence c = (ci, . . . , Cm) G ffi > Ci > ■ ■ ■ > Cm > 0, denote by 
f2^(-F) C Gr(y, M) the corresponding Schubert cell. 

For i — 1, . . . , N, define the ramification sequence 

c(i) = (rii, 0, . . . , 0) . 

Define the ramification sequences 

c(cx)i) = (m2 H h rriM, ^3 H h uim, ruM, 0) , 

0(002) = {mi-\ \-mM-i:mi-\ h mM-2, • • • , ""^i, 0) . 

Lemma 7.6. 
• We have 

N 

J2 codimn^(,)(F(A,)) + codimn:;(^^)(F(ooO) + 
1=1 

codim ^(F(oo2)) =dim Gr(y, M) = Mm . 
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• Let U E Gt(Y,M). The pair {U, X) is a space of the {m, z,n, X)-type, if and 
only if U belongs to the intersection of N + 2 Schubert cells 

fi:(i)(F(Ai))nr]°(2)(i^(A2))n... 

□ 

According to Schubert calculus, the multiplicity of the intersection of Schubert cycles 

f]e(l)(i^(Al))nfie(2)(i^(A2))n... 

(7.4) n n,(^M)iF{\N)) n n^^^^){F{ooi)) n n,(^^^){F{oo2)) . 

can be expressed in representation-theoretic terms as follows. 

For a ramification sequence c = {ci, . . . ,cm) E Z*^, denote by Li^^ the finite di- 
mensional irreducible gl^^-module with highest weight c. Any glj^^-module L^*^^ has a 
natural structure of an slAf-module denoted by L^^'^K By |Fuj . the multiplicity of the 
intersection in 1)7.11) is equal to the multiplicity of the trivial stM-module in the tensor 
product of slM-niodules 

By Proposition [TiSl the multiplicity of the trivial slM-module in the tensor product ()7.5p 
is equal to the dimension of the weight subspace of weight [mi, . . . , mjv/] in the tensor 
product of gljvf -modules 

(7.6) 

Corollary 7.7. For generic z, the number of orbits of the XXX admissible of X-type 
solutions of the Bethe ansatz equations \6.S\) is not greater than the dimension of the 



weight space ® • • • ® ^^c(n))[^^i • • • ' "^m] 



□ 



8. The KZ and dynamical Hamiltonians 



8.1. The Gaudin KZ Hamiltonians. Let Eij, i,j = 1,. . . ,N, be the standard gen- 
erators of the complex Lie algebra gl^. 

We have the root decomposition gt^ = n"*" © f) © where 



)i<j C ■ Eij , () — ©i=l ^ ' -^ii 1 ^ — ®i>j C ■ Eij 



Let 

N 



= ^ ^ ® Eii , = Q^ + Eij © Eji , Q- = + Eji © Eij . 

i=l i<j i<j 

Let Y = Fi©. . .©Ym be the tensor product of finite-dimensional irreducible gt^-modules. 
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The Gaudin KZ Hamiltonians Hf{\, z), a = 1, . . . , M, acting on y-valued functions 
of A = (Ai, . . . , Atv), z = {zi, . . . , zm) are defined by the formula |TV5j : 

■ 1 ^ I, 1 Za Zl) 

1=1 0=1, b^a 

Here the hnear operator -.Y^Y acts as on Ya Vf,, and as the identity on 

other tensor factors of Y . Similarly, E^f^ acts as En on and as the identity on other 
factors. 

8.2. The Gaudin dynamical Hamiltonians. For any i, j = 1, . . . , N, i j, intro- 
duce a series Bij{t) depending on a complex number t\ 

oo s ^ 

The series has a well-defined action in any finite-dimensional gt^-module W giving an 
End (Vr)-valued rational function of t. 

The Gaudin dynamical Hamiltonians Gf{X,z), i = 1,...,N, acting on F- valued 
functions of A, z are defined by the formula |TV5j : 

Gf{X, z) = {Bi^j^{Xi — Xn) ■ ■ ■ -Bi,i+i(Ai — Aj+i))"^ x 

M 

H X i?i,(Ai - A,) . . . i?.-i,.(A,_i - A,) . 

a=l 

8.3. The Gaudin diagonalization problem. The Gaudin KZ and dynamical Hamil- 
tonians commute ^T V4^ . 

[H^iX,z),H^iX,z)] = 0, [Hf{X,z),Gf{X,z)] = 0, [GfiX,z),Gf{X,z)] = 0, 

for a, 6 = 1, ... , M, and i,j = l,...,N. 

The Gaudin diagonalization problem is to diagonalize simultaneously the Gaudin KZ 
Hamiltonians Hf, a = 1, . . . , M, and dynamical Hamiltonians Gf, i = 1, . . . ,N, for 
given A, z. The Hamiltonians preserve the weight decomposition of Y and the diagonal- 
ization problem can be considered on a given weight subspace of Y. 

8.4. Diagonalization and critical points. For a nonnegative integer m, denote by 
Lm'^ the irreducible finite-dimensional gl^y-module with highest weight (m, 0, . . . , 0). 

Let n = (ni, . . . , n^) and m = (mi, . . . , rriM) be vectors of nonnegative integers with 

EN v-^A/ 

Consider the tensor product Lm^ ® . . . Lm}j and its weight subspace 
Lm[n] = (lW®...®L<„^2)K'---'^^]- 
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Let A G C^, z G C^^. Assume that each of A and z has distinct coordinates. Consider 
the Hamiltonians Hf{X,z), Gf{X,z) acting on L^[n]. The Bethe ansatz method is a 
method to construct common eigenvectors of the Hamiltonians. 

As in Section 16.11 consider the space C"^"* i-"]v-i -with coordinates t^"'^ Let 
A; m; z) be the master function on +"^-1 defined in ()6.1|) . 

In Section 4 of |MaVj . a certain L^[n]-valued rational function tu® : ([y^+-+'<^N^i _^ 
Lm[ri\^ depending on z, is constructed. It is called the Gaudin universal rational func- 
tion. For N = 2 formulas for the Gaudin universal rational function see below in Section 

ino 

The Gaudin universal rational function is well defined for admissible t^""^ . The Gaudin 
universal rational function is symmetric with respect to the Xlf^-action. 

Theorem 8.1 ( |RV[ IMaVj ). //t^"^ is a Gaudin admissible non- degenerate critical point 
of the master function $(• ; A — n; m; z), then Lu'^{t^'^\ z) G -LmM is an eigenvector of 
the Gaudin Hamiltonians, 

(8.1) iff(A,z) c.*(t<-\z) = ^ log A -n;m;z) a; V\ ^) , 

a = 1,...,M , 

M ni 

a=l j=l 
fii-i fii 

Gf{\z)u^{t^-\z) = l[(t^-'Y'lltfu;^it^,z), ^ = 2,...,iV-l. 

njv-i 

G^(A,z) co^{t^-\z) = n i^r'^r' u^^it^'-Kz) . 

If t^"'^ is a Gaudin admissible critical point of the master function, then the vector 
uj^(t^'^\z) is called a Gaudin Bethe eigenvector. 

By Corollary 17. 4t for generic A, the number of Gaudin Bethe eigenvectors is not 
greater than the dimension of the space 2y^[n]. The Bethe ansatz conjecture says that 
all eigenvectors of the Gaudin Hamiltonians are the Gaudin Bethe eigenvectors in an 
appropriate sense, cf. |M V2j . 

8.5. The XXX KZ Hamiltonians. Let Eat, a, 6 = 1, ... , M, be the standard gener- 
ators of the complex Lie algebra qIj^j. 

We have the root decomposition 01^ = ^^®') © ri^ where 



n+ = ®a<b C-Eab, f) = ©^ii C-Eaa, U" = (Ba>b C ■ E, 



ab ■ 



Let V,W be irreducible finite-dimensional glj^^-modules with highest weight vectors 
V E V, w E W. The associated rational R-matrix is a rational End (V ® Vr)-valued 
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function Rywit) of a complex variable t uniquely determined by the g[^-invariance 
condition, 

[Rvw{t) , 9(^l + l(» g] = for any g e , 
the commutation relations 

M M 

Rvw{t) {tEal,®l + ^Eac<^ E^b) = {tEij®l + J2^cb®Eac) Rvw{t) , 

c=l c=l 

for any a, b, and the normalization condition 

Rvwit) V 1^ w = V0W. 

Let Y = Yi ^ . . . Yjy he the tensor product of finite-dimensional irreducible gl^j- 
modules. 

The XXX KZ Hamiltonians Hf{z, X), i = 1, . . . , N, acting on F- valued functions of 
z = {zi, . . . , zm), a = (Ai, . . . , Xn) are defined by the formula |TV5j : 

Hf{z, A) = (-Ri,Ar(Aj — Aat) . . . Ri^i+i{Xi — Aj+i))"^ x 

M 



H (^-^"»)(*) X i?i,,(Ai - A,) . . .i?,_i,,(A,_i - A,) 



a=l 



where Rjk{t) is the endomorphisms of Y which acts as Ry^Vkit) ® Y^ and as the 

identity on other factors of Y. 

8.6. The XXX dynamical Hamiltonians. The XXX dynamical Hamiltonians 
G^{z,X), a = 1, . . . , M, acting on F-valued functions of z,X are defined by the for- 
mula [TV5] : 



,2 ^ 



2 

i=l 



M M 

J2 -^{EabEta-Eaa) + (i^afc) ^^^i^ba) ^^'^ 

b=l, b^a b=l i<j 



8.7. The XXX diagonalization problem. The XXX KZ and dynamical Hamiltoni- 
ans commute |T V6j . 

[H^{z,X),Hf{z,X)] = 0, [H^{z,X),G^{z,X)] = 0, [G^{z, X),G§{z, X)] = 0, 

for i,j = l,...,N, and a, 6 = 1, ... , M. 

The XXX diagonalization problem is to diagonalize simultaneously the XXX Hamil- 
tonians for given z, X. The Hamiltonians preserve the weight decomposition of Y and 
the diagonalization problem can be considered on a given weight subspace of Y. 
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8.8. Diagonalization and solutions of the Bethe ansatz equations. For a non- 
negative integer n, denote by lI^^'* the irreducible g [^-module with highest weight 
(n,0,...,0). 

Let m = {nil, ■ ■ ■ , ^im) and n = {ni, . . . , un) be vectors of nonnegative integers with 

EM sr^N 

Consider the tensor product lI*^^ ... lI^^ and its weight subspace 



Let z G C^, A G C^. Assume that each of z and A has distinct coordinates. Consider 
the Hamiltonians H^{z,X), G^{z,X) acting on Ln[m]. The Bethe ansatz method is a 
method to construct common eigenvectors of the Hamiltonians. 

As in Section ini21 consider the space C^i+ --+'"m-i -^^ith coordinates t^"^^ and the XXX 
Bethe ansatz equations (|6.3|) . Let 1 = (1, . . . , 1) and 



\^-^) S ISl ) ■ • ■ ) Smi ) SI )---)Sm2? SI ) • • • ) SmAf_i y 



be the point with coordinates ^ = for i = 1, . . . , M — 1 and a = 1, . . . , rfiM-i- 
In |TVlj . a certain Zy^ifm] -valued rational function : C"^^"^ ""*"'"""^ ^ -Lri[m] , 

depending on A , is constructed. It is called the XXX universal rational function. For 

N = 2 formulas for the XXX universal rational function see below in Section [10.31 
The XXX universal rational function A) is well defined if is XXX 

admissible of {X — n + l)-type. The XXX universal rational function is symmetric with 

respect to the S^-action. 



Theorem 8.2. [TV2llMTV3l //t<"^)+^("^> is an XXX admissible solution of{X-n + l)- 
type of the Bethe ansatz equations \6. 'J\} . then uj m] is an eigenvector of 

the XXX Hamiltonians, 



mi ^(1) _ ^ 

(8.3) H^{z,X)ujHt^^\X) = HlTr —uj^t^'^KX), z = l,...,N, 

fj^ tf - Xi + n. 



Gf{z, A) u;^{t^^\ A) = d,{t^^\X) u;^{t^^\ A) , a = 1, . . . , M , 
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where 

N ?fii _ M 





i=l 


1 + E 4" 


mi 
2" 


-E 

6=2 


rribZb 

Zi — Zb 


fria 


- E 'r" - 




a-l 

-E 

6=1 


rriaZb 


M 

- E 

6=a+l 


2 





rribZb 

Za - Zb 

a = 2,...,M-l 



niM-i _ A/-1 

(a-l) rriM-i "^4/^6 



<iM(t<"->,A) = - ^ 4°-" - ^ - E 



j=i 6=1 

If + is an XXX admissible solution of (A - n + l)-type of the XXX Bethe 
ansatz equations, then the vector uj^{t^'^\ A) is called an XXX Bethe eigenvector. 

By Corollary 17.71 for generic z the number of the XXX Bethe eigenvectors is not 
greater than the dimension of the space 2y„[m,]. The Bethe ansatz conjecture says that 
all eigenvectors of the XXX Hamiltonians are the XXX Bethe eigenvectors. 

According to previous discussions, eigenvectors of Gaudin Hamiltonians on -Lrn[*^] are 
related to critical points of the master function A — n; m; z). The critical points 

of A — n; m; z) are related to spaces of the (n, A — n, m, z)-type. The spaces of 

the (n. A — n, m, z)-type are special bispectral dual to spaces of the (m, z, n. A — n + 1)- 
type, where 1 = (1, . . . , 1) G . The spaces of the (m, 2, n, A — n + l)-type are 
related to solutions of the XXX Bethe ansatz equations, which in its turn are related 
to eigenvectors of the XXX Hamiltonians acting on Iv„[m]. As a result of this chain of 
relations, the eigenvectors of the Gaudin Hamiltonians on and the eigenvectors 

of the XXX Hamiltonians on X^ii^iT'] must be related. 

Indeed this relation is given by the {qI^ , duality. 

9. The (gl^ .Qhi) duality for KZ and dynamical Hamiltonians 

9.1. The (g[^ , qIm) duality. The Lie algebra qI^ acts on C[a;i, . . . , xat] by differential 

d 

operators E-ij ^ Xi— — . Denote this gl^y'^odule by L,. Then 



00 

{N) 



m ' 

m=0 



the submodule Lm'^ being spanned by homogeneous polynomials of degree m. The gi 
module Lm'^ is irreducible with highest weight (m, 0, . . . , 0) and highest weight vector 

-^1 • 

We consider gl^ and gl^^ simultaneously. To distinguish generators, modules, etc., we 
indicated the dependence on and M explicitly, for example, , Lm^\ 



■N' 
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Let Pmn = C[xii, . . . , xmi, • • • , xin, • • • , xmn] be the space of polynomials of MN 
variables. We define the g [^-action on Pmn by -E^^^^ ^ — gfTV-action 

i=l 

by E^^'^ t-^ — ■ There are two isomorphisms of vector spaces, 

a=l ■' 

(9.1) 

N 

(C[xi, . . . ,xm])^^-^ Pmn, (pi ® . . . ® PAf)(a;ii, . . .,xmn) ^ W Piixu, . . .,XMi), 

i=l 

M 



Xi, 



xat])® ^ Pmn, (pi ® • • • ® PM)(a;ii, . . . , xma?) J]^ Pa(a;ai, • • • , a;aAr). 



a=l 



Under these isomorphisms, Pmn is isomorphic to (^{^[^^'^^^^ as a glj^^-module and to 
(L.<^>)®^ as a sljv-module. 

Fix n = (ni, . . . , n^) G ^nd m = (mi, . . . , ttia/) G Z^^, with Xl^i «i = Zlali "^a- 
Isomorphisms ()9.ip induce an isomorphism of the weight subspaces, 

(9.2) L^[m] ~ X^[n] . 

Under this isomorphism the KZ and dynamical Hamiltonians interchange, 
Hf'(^\X,z) = G^'^'^\z,X), G?'<^^(A,z) = //f'<^H-,A), 
for a = 1,...,M, i = l,...,iV, lTV5] . 



Recall that 1 = (1, . . . , 1) G C^, and is defined by 
Let t^"^ be a Gaudin admissible critical point of the master function 

$( ■ ; A — n; m; 5;). Let (V, z) be the associated space of the (n, A — n, m, 2)-type. Let 

{U, X — n + 1) be its special bispectral dual space of the (m, z,n, X — n + l)-type. 

Assume that (f/, A - n + 1) is XXX admissible. Let + ^^"'^ be the associated XXX 

admissible solution of the Bethe ansatz equations ()6.3|) . 

Conjecture 9.1. The corresponding Bethe vectors uj'^{t^'^\z) G Lm[Ti] and 
uj^(t^'^\ X) E Ln[m] are proportional under the duality isomorphism 119.^) . □ 



10. The case of = M = 2 

10.1. The {q12,Q^2) duality. Let n = (ni,n2) and m = (mi, 7/12) be two vectors of 
nonnegative integers such that ni + n2 = mi + m2. 

Fix A = (Ai, A2) and z = (zi, Z2) each with distinct coordinates. 

For m G Z>o, let Lm be the irreducible gl2-in^odu\e with highest weight (m, 0) and 
highest weight vector Vm- The vectors E'aifm, i = 0, . . . , m, form a basis in L^- 



BISPECTRAL AND (bIjv.sIm) DUALITIES 39 

Set 

a = max (0, n2 — rrii) , j3 = min (m2, . 

The vectors 

form a basis in Lmln] = {Lmi ® Lm2)[ni,n2]. The vectors 

7 TT® — ^ — , a<«<P, 

form a basis in 'LrSrn\ = {Ln^ (g) L„2)[mi, 7712]. Isomorphism ()9.2|1 identifies the vectors 
with the same index i. 

We consider the Gaudin Hamiltonians on -Lrn,[^] and the XXX Hamiltonians on 
Ln[m\. 

Let Bij{t) be the series defined in Section IH^ Consider the map 

(10.1) 6(A) : (L„^ (g) L„J[ni,?22] -> (L^^ ® L„J[n2, ^1] 
defined by the formula 

7 ^ ® 1 — ^2UAi - As) —7 ® 7 — - . 

For a generic A this map is an isomorphism. 

By |TV4j . the map 6(A) commutes with the Gaudin Hamiltonians in the following 
sense: 

b{X)H^iXuX2,z) = i7f(A2,Ai,z)6(A) , 
6(A)Gf(Ai,A2,z) = Gf(A2,Ai,2)6(A) 

for a,i = 1, 2. 
The map 

J?m2-i pi pni-m2+i R"2-i 

(10.2) y ^n. ^^^^ ^ ^^f:2L^ 

(m2 — i)! V. (ni — m2 + i)' \n2 — '^)- 

defines the Weyl isomorphism 

(10.3) s : (L„j (g) L„2)[mi,m2] (L„^ (g) L„J[m2, mi] . 
The Weyl isomorphism commutes with the XXX Hamiltonians, 

sEf{z,\) = H^{z,X)s, 

for i,a = 1, 2. 
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10.2. The Gaudin universal rational function and the fundamental differential 
equation. Let n = {ni,n2) and m = (1711,7712) be two vectors of nonnegative integers 
such that 721+^2 = mi + 

Fix A = (Ai, A2) and z = [zi, Z2) each with distinct coordinates. The Gaudin universal 
rational Ljn[n]-valued function is the function 

Tpn2—i rpi 
1, J- ^ _ \^ n 21 '^mi ^ ^2\^m2 

U)\t\, . . . , tn^) — 



[n2-ty. %\ 



where 



n2-t - ^ 

(10.4) Q(ti, . . . , t„J = Sym„^ J] 7^ H 7 

and Sym„/(ti, . . . = Eags„ • • • ,^a(n)), see |MaVj. 

The Gaudin universal rational function is symmetric with respect to the group 
of permutations of variables ti, . . . ,t„2- The S„2-orbit of a point t^"'^ = (ti, . . . ,)f:„2) is 
represented by the polynomial p{x) = {x — ti) . . . {x — tn^) . 

The polynomials {x — zi)"^~'^{x — ^2)*, = 0, . . . ,n2, form a basis in the space of 
polynomials in x of degree not greater than n2. 

Lemma 10.1. Let p{x) = (x — ti) . . . (a; — t„2) ^ polynomial. Let the numbers 
Co, ... , given by formula lilU.4\ )- Then 

p{zi)p{z2) ^ * (n2-«)! i\ 

□ 

Let V = {x^'^pi{x) , x^^p2{x)). Let (V^, z) be a space of the (n, A, m, z)-type. The 
special fundamental differential operator of {V, z) has the form 

D = {x ~ zi){x - Z2){xdr, - Xi - ni){xdj: - X2) 

+ x{x - zi){xd^ - Ai - ni) + 012 Z2 {x - zi){xd^ - A2) 

+ 021 x{x - Z2){xdx - Ai - rii) + 022 zi {x - 2:2) - A2) , 

where 0jj are suitable numbers such that 

(10.5) 021 + 022 = -mi, 011+012= -m2, 

022 + 012 = - rii, 021 + 011 = - n2. 

The operator D can be written also in the form 

D = {x - zi){x - Z2){xd^ - Xi - ni){xd^ - X2 - 712) 

+ ^11 X{X - Zi){xd:r " Ai) + V^12 Z2 {x - Zi){xdx - A2 - 712) 

+ 1^21 x{x - Z2){xdx - Ai) + V^22 Zi {x - Z2){xdx - A2 - n2) , 
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where ijjij are such that 



4^21 + ip22 = -mi, 
^22 +^12 = -ni, 



^21 + ^11 = -n2, 



and ipu{Xi — X2 



"-2) = 01l(Al - A2 + rii) + 7712^2 . 



Equation Df 



has a solution 



a; ^^2(2;) 



J2 



(x - 2:1) 



= X' 



i=0 



n2 - t)l 



Due to conditions (jlO.Sp . the expression x^^^Df{x) is a polynomial of degree at most 
77,2 + 1 vanishing at x = 0. Expanding this polynomial as a linear combination of the 
polynomials x{x — zi)^^^'^{x — 22)% i = 0, . . . , n2, we obtain that the equation Df = is 
equivalent to the following relations for the coefficients cq, . . . , Cnj: 

(10.6) zii{mi - n2 + i) Ci-i + Z2{m2 - i){n2 - i) Ci+i 

+ ((2:1 + Z2) + {z2mi - Zim2 - 2z2^2 - (Ai - A2)(2;i - Z2)) i 

- (Ai - A2 + ni){zi - 2:2)011 + Z2m2n2) Q = , 

i = 0, . . . ,n2. Notice that values i = 0,722 — mi, 7712,722 are the values of i for which 
equation p0.6|) does not contain q_i or q+i. 

Lemma 10.2. Equations ilU.6]) with i such that a < i < j3 form a dosed system of 
equations with respect to cj such that a < j < (3. □ 

Equations p0.6p have a symmetry. Namely, ()10.6p does not change if we replace the 
parameters Ai, A2, tii, 712, 0ii, ^ by A2, Ai, 7Z2, 77-1, ?/'i2, m2 — ^, respectively, and after that 
replace the unknowns Cj by {zi/ Z2y Cm2-j- 

10.3. The XXX universal rational function and the fundamental difference 
equation. Let n = (711,772) and m = (7771,7772) be two vectors of nonnegative integers 

such that 77i + 772 = 77li + 77^2. 

Fix A = (Ai,A2) and z = {zi,Z2) each with distinct coordinates. The XXX universal 
rational Ln[m]-valued function is the function 




42 

where 

(10.7) a(si,...,s^„A) 
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Sym 



m2 



1712-1 



n 



\ Sa- \i + ni 



m2 

n 

b=m2—i+l 



Sb - Ai 

(sb - Ai + ni) {sb - A2 + 



X 



m2—i 1712 

n n 



- - 1 

- Sb 



a=l 6=m2— 2+1 

see |KBTHTV2] . 

The XXX universal rational function is symmetric with respect to the group of 
permutations of variables Si, . . . , s^j. The Smj-orbit of a point s^"'^ = (si, . . . , s^a) is 
represented by the polynomial q{u) = {u — Si) . . . {u — s^j). 

The polynomials 11^^(7*"^ - Ai + ni - j) nj=o(^ - A2 + j) 5 = • • • > "^2, form a 
basis in the space of polynomials in x of degree not greater than 1712. 

Lemma 10.3. Let q{u) = {u — Si) . . . {u — s^j) be a polynomial. Let the numbers 



Co, 



, Cm2 be given by formula ^lU. 7\ ). Then 

7712 — 1 

JJ (Ai - A2 + n2 - j] 



q{u] 



g(Ai - ni) g(A2 - 712) ^ 



j=0 



Si, . . . ,s 



) ''ma ) 



A) 



m2— 2— 1 



j-1 



j=0 



(m2 — tjltl 



Yl (n - Ai + j) JJ (n - A2 + ^2 - j] 



j=0 



j=0 



□ 



Let U = Z2q2{u)). Let (f/, A) be a space of the (m, 2, n, A)-type. The 

special fundamental difference operator of (t/, A) has the form 

D = {u- Xi-ni + l){u- X2 + 1){t - zi){t - Z2) 

+ (pu {u - Xi - rii + l)r(r - Zi) + ipi2 {u - X2 + 1)-22(t" - Zi) + 
+ '-P21 (u - Ai - + l)r(r - 2:2) + ip22 (m - A2 + l)2;i(r - Z2) , 
where (/^jj are suitable numbers such that 

(10.8) ^P2i + '■P22 = -mi, V^ii + = - ^2, 

<^22 + <-Pi2= - ni, ip21 + <-Pii= - n2. 

Equation Df = has a solution 

1712 1712—1 — 1 i—1 

f = zl^q2{u) = z- n («-Ai-ni+j) l[(u-X2-j) . 



m2 — 2 u 

i=l ^ ^ ^ j=0 
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Due to conditions p().8|l . the expression Z2^Df{u) is a polynomial of degree at most 
m2- Expanding this polynomial as a linear combination of the polynomials 

m2— j— 1 i—1 

Y\ (m - Ai - ni + j + 1) J]^ (m - A2 - j + 1) , i = 0, . . . ,1712 , 

we obtain that the equation Df = is equivalent to the following relations for the 
coefficients cq, . . . , Cm2- 

(10.9) Zii{mi - n2 + i) Ci-i + Z2{m2 - i){n2 - i) Ci+i 

+ [{zi + Z2) + {z2mi - Zim2 - 2z2n2 - (Ai - \2){zi - Z2)) i 

- (Ai - A2 + ni){zi - Z2)(pii + Z2m2n2) q = , 

i = 0, . . . , m2. Notice that values i = 0,n2 — 1711,1712, 712 are the values of i for which 
equation pO.fjj) does not contain q_i or 

Lemma 10.4. Equations MO.f^) with i such that a < i < j3 for7n a closed system of 
equations with respect to cj such that a < j < jS . □ 

Equations ()10.9p have a symmetry. Namely, ()10.9p does not change if we replace the 
parameters zi, Z2,mi,m2, (pn, i by ^2,-21, ^2, fni, ip2i,n2 — h respectively, and after that 
replace the unknowns cj by c„2_j. Another symmetry of ()10.9|1 is that the equations do 
not change if we replace Ai, A2 by Ai + 1, A2 + 1. 

The pair of equations (jlU.6p and (jl(J.9j) have a symmetry. Namely, equation (|1(J.6|) 
turns into equation (jl(J.9|) if we replace the parameter (pu there by ipu. 

10.4. Proof of Conjecture 19.11 for N = M = 2. The symmetries of equations (jlO.fij] 
and (jl0.9p imply the following theorem. 

Theorem 10.5. Let V = {x^^pi{x) , x^'^'^p2{x)) and U = {ziqi{u) , z'2q2{u)) . Let 
{y, z) be a space of the (n. A, m, z)-type and {U, A + 1) the special bispectral dual space 
of the (m, z,n, X + l)-type, where 1 = (1, 1) . Write 



n.2 



P2[X) 



X — * (x — Z2) 



ni 



nil f ( \ _i_ -W mi-i-l i-1 

*w = E7^riyT n (''-^■-i-".+j-)ne"-^=-i-j) 

i=l ^ ^ ^' ' j=0 j=0 
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Then there exist nonzero numbers A{X), B{X),C{X) such that 

c,(A) = A{X)dm,-^{X){zi/z2y = B{X)e,{X + l) = C(A) /„,_,(A + 1) 
for all i, a < i < (3. 

Theorems 110.51 and 18. II yield the next statement. 
Corollary 10.6. The Bethe vectors 

^-^ (mo — ty. il ^-^ (mi — zi! z! 

are identified up to proportionality by isomorphisms ilU.l}} and OJ.!^) . □ 

The fact that the first and second Bethe vectors are identified by isomorphism po.lj) 
is proved also in ^MV3) in a different way. 

10.5. Correspondence of solutions of the Bethe ansatz equations for N = M = 

2. Consider the following four systems of the Bethe ansatz equations: 

(10.10) / + Er^V - E = 0' ^ = 1' 

a=l j=l, j^i -J 



,n2 



2 _ ni 



'10.11) ^1^^ 



a=l j=l, j^i 

~ Aj — 1 T-r ta ~ tf) — 1 Z2 



'10.12) n " n 

' J-J- - A, - l-n,. IJ- 



m2 



■ 1 - - ^ I, — h + ^ Zi 

1=1 0, b=ia 

(10.13) n ^--^--^ fj ^"-^^-1 = a = l,...,mi. 

The first two systems are of the Gaudin type and the last two systems are of the XXX 
type, see (Q and (jHISI)- 

Set d = dim{Lrn^®Lrn2)[ni,n£= dim{Lrn^®Lrn2)[n2,ni] = dim (L„^(g)L„J[mi, = 
dim (L„^ (g) L„2)[m2,mi] . By )MV4j for generic A and z each of the first two systems 
has exactly d orbits of Gaudin admissible solutions. By |MV4j for generic A and z each 
of the last two systems has exactly d orbits of XXX admissible solutions. 

Theorem 10.7. Let X and z be generic. Let {ti, . . . , tnj) be a Gaudin admissible solution 
of the system ilU.lU^) . Let (V, z) = {x^^pi{x),x^'^p2{x)) be the space of the {X, z,n,m)- 
type associated with the orbit of the critical point. Let {U,X+ 1) = {z^qi{u), Z2q2iu)) , 
where 1 = (1,1), be the special bispectral dual space to (y,z). Then the polynomial pi 
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represents the orbit of a Gaudin admissible solution of the system M().ll\) . the polynomial 
q2 represents the orbit of an XXX admissible solution of the system ilU.l^) . and the 
polynomial qi represents the orbit of an XXX admissible solution of the system MU.1'J\} . 

□ 



11. BAKER-AkHIESER functions AND BISPECTRAL CORRESPONDENCE 

11.1. Grassmannian of Gaudin admissible non-degenerate spaces. For A G C, 

we call a vector subspace Xx C C[m] Gaudin admissible at A if there exists m G Z>o 
such that 

(u - A)(m - A - 1) . . . (m - A - m + 1) C[u] C Xx 

and there exists / G Xx such that /(A) ^ 0. We call a vector subspace X C C[?i] Gaudin 
admissible if 



A- = n X- 



1 



where n is a natural number, Ai, . . . , A„ G C are such that Aj — Aj ^ Z if i 7^ j, and for 
any i, Xx^ is Gaudin admissible at Aj. 

Each Xxi can be defined by a finite set of linear equations. Namely, there exist a 
positive integer Ni, a sequence of integers Q < rin < ■ ■ ■ < UiN^ and complex numbers 
Ci,j,a5 a = 0, . . . , riij, such that Cjj^„.^ 7^ for all z, j. Then the space X^. consists of all 
polynomials r G satisfying the equations 



;il.l) ^ c,,,- , r(A, + a) = for j = 1, . . . , iV, 



a=0 

For a Gaudin admissible subspace X, define the complex vector space V as the space 
spanned by functions x^^pij{x), i = 1, . . . ,n, j = 1, . . . , Ni, where Pij{x) = J2a=o Cij^aX""- 
The space is a space of quasi-polynomials of the type considered in Section ITTl 

On the other hand, having a space V of quasi-polynomials, as in Section 12. H we can 
recover a Gaudin admissible subspace X C C[u] by formula (jll.ip . 

We say that the space X is non- degenerate if the corresponding space V is non- 
degenerate in the sense of Section 12.11 

We denote the set of all Gaudin admissible non-degenerate subspaces by Gr® and call 
it the Grassmannian of the Gaudin admissible non-degenerate subspaces. 

Let X C C[n] be a Gaudin admissible subspace and V the associated space of quasi- 
polynomials in X. Define the algebra Ax = {p G <C[u] \ p{u) X G X}. An equivalent 
definition is Ax = {p G C[u] \ p{xdx) V dV }. 

Let Dy be the monic fundamental differential operator of V . The function 

n Ni 

= nn 1731^ 5,." 
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is called the stationary Baker-Akhieser function of the Gaudin admissible space X. In- 
troduce the rational function ■iIjx{x,u) by the formula '^x{x,u) = ipx^x^u) . The 
function ipxi^iu) expands as a power series in of the form 

oo 

1px{x,u) = 1 + ^ CijX"' U"^ . 

It is easy to see that for every p G Ax, there exists a linear differential operator Lp{x, dx) 
with rational coefficients such that Lp{x,dx) Dy = Dvp{xdx). As a corollary, we con- 
clude that 

Lp{x,dx)^!x{x,u) = p{u)'^x{x,u) . 
For pi,p2 G A, the corresponding operators Lpj^{x,dx) and Lp^{x,dx) commute. 

11.2. Grassmannian of XXX admissible non-degenerate spaces. Let 2 G C* be 

a nonzero complex number with fixed argument. We call a vector subspace C C[x] 
XXX admissible at z if there exists n G Z>o such that {x — z)"^ C[a;] C Yz and there exists 
f &Yz such that f{z) ^ 0. We call a vector subspace F C C[x] XXX admissible if 

m 

^' = n y^. 

a=l 

where m is a natural number, zi, . . . , G C* are distinct numbers and for any a, the 
space Yz^ is XXX admissible at Za. 

Each Yz^ can be defined by a finite set of linear equations. Namely, there exist a 
positive integer M^, a sequence of integers < niai < ■ ■ ■ < rriaMa complex numbers 
da,b,i, i = 0, . . . , niab, such that da^b,mab 7^ ^- Then the space consists of all 

polynomials r G C[a:;] satisfying the equations 

mat 

(11.2) 5^ cia,M^^'^(^a) = for b = l,...,Ma. 

i=0 

For an XXX admissible subspace Y define the complex vector space U as the space 
spanned by functions Qatiu), a = 1, . . . , m, 6 = 1, . . . , Ma, where 

rriab 

^aQabiu) = ^ ^ da,b,i X \x=Za ■ 

To determine x" at x = Za we use the chosen argument of Za- 

The space f/ is a space of quasi-exponentials in u of the type considered in Section 

o 

On the other hand, having a space U of quasi-exponentials, as in Section 13.11 we can 
recover an XXX admissible subspace F C C[x] by formula ()11.2|) . 

We say that Y is non- degenerate if the corresponding U is non-degenerate in the sense 
of Section 13.61 
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We denote the set of all XXX admissible non-degenerate subspaces by Gr^ and call 
it the Grassmannian of the XXX admissible non-degenerate subspaces. 

Let Y C C[a;] be an XXX admissible subspace and U the associated space of quasi- 
exponentials. Define the algebra Ay = {p G C[a;] | p{x) W C W}. An equivalent 
definition is Ay = {p G C[x] | p(r„) U C U}. 

Let Djj be the monic fundamental difference operator of U . The function 

m 
a=l 

is called the stationary Baker- Akhieser function of the XXX admissible space Y. Intro- 
duce the rational function 0y(u, x) by the formula x) = (f)y{u,x) x"^. The function 
4>y{u, x) expands as a power series in u~^, x~^ of the form 

oo 

(j)y{u,x) = 1+ CijU~^X~^ . 

It is easy to see that for every p G Ay, there exists a linear difference operator Lp{u, r^) 
with rational coefficients such that Lp{u, r^) Dy = Dvp{tu). As a corollary, we conclude 
that 

Lp{u,ru) ^y{u,x) = p{x)^y{u,x) . 
For pi,p2 G A, the corresponding operators Lp^{u,Tu) and Lp^{u,Tu) commute. 

11.3. Bispectral correspondence. In terms of the Baker-Akhiezer functions, we in- 
troduce a bispectral correspondence between points of the Grassmanian of Gaudin ad- 
missible non-degenerate subspaces and points of the Grassmanian of XXX admissible 
non-degenerate subspaces. We say that a Gaudin admissible non-degenerate space X 
corresponds to an XXX admissible non-degenerate space Y if "^xix, u) = $y(m, x). This 
correspondence is an analog of Wilson's bispectral correspondence in jWj. 

Theorem 11.1. Let V be a non-degenerate space of quasi-polynomials and let U be the 
non- degenerate space of quasi- exponentials which are bispectral dual with respect to the 
integral transforms of Section [T} Then the corresponding admissible spaces X G Gr*^ 
and Y G Gr^ are bispectral correspondent. □ 
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